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PREFACE 


This  study  vas  doxw  at  the  Visibility  Laboratoxyi  Soripps 
Institution  of  Oceanography^  University  of  California^  San 
Diego,  by  Mr.  Richardson  outside  of  vorking  hours  as  his  thesis 
for  the  degree  of  Master  of  Science  in  Engineering  from  the 
Urdverslty  of  California,  Los  Angeles.  Although  the  study 
was  done  at  no  cost  to  the  Laboratory  oontracts,  it  is 
pertinent  and  of  interest  to  the  program  of  the  Laboratory 
and,  as  such,  is  presented  as  a  Laboratory  report. 


S.  Q.  Duntley 
Director 

Visibility  Laboratory 
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ABSTRACT  OF  THE  THESIS 


A  Study-  of  tho  Factors  Affaoting  the  Sighting 
of  Surface  Vesaela  from  Aircraft 

by 

William  Hadley  Rlohardaon 

A  oolleotlon  of  3|465  detailed  reports  of  alghtlnga  of  surface 
vesaels  from  aircraft  of  the  U.  S.  Coast  Guard  are  analyzed  by  pro¬ 
bit  analysis  to  determine  visual  thresholds  and  measures  of  -varianos 
of  the  thresholds.  Each  of  seventeen  conditions  affecting  the  sight¬ 
ing  range  is  studied  separately  to  determine  its  effect.  Easpirioal 
funotions  are  developed  to  describe  the  threshold  effects  of  each  of 
the  eleven  following  conditions  in  decreasing  order  of  inportancei 
(1)  meteorological  visibility,  (2)  altitude  of  aircraft,  (3)  ship 
sise,  (4)  height  of  major  swells,  (^)  cloud  cover,  (6)  wind  velocity, 
(7)  relative  bearing  of  teurget,  (8)  sun  altitude,  (9)  relative  bear¬ 
ing  of  a-on,  (10)  wake  size,  and  (11)  wind  azlmu-th.  Thresholds  are 
developed  for  the  six  followitig  discrete  conditions  in  decreasing 
order  of  importance:  (12)  visual  aid,  (13)  range  determination 
method,  (14)  type  of  observing  unit,  (15)  time  of  day,  (16)  obsetver, 
(17)  station,  in  each  case,  tables  of  problt  results  and  graphs  ars 
Included,  Only  wind  azimuth  is  found  to  have  an  insignifioant  affect. 
Classification  of  data  is  made  by  mechanical  card  sorter,  probit 
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analysis  by  electronic  computer  and  the  remainder  of  the  calcula¬ 
tions  by  desk  oaloulator.  Measures  of  precision  are  listed  and 
X  ,  F  and  t  tests  are  used  at  the  0.05  level.  A  table  of  fac¬ 
tors  for  each  condition  is  included  to  allow  forecasting  of  sighting 
thresholds  and  explanation  of  use  for  any  probability  level.  A 
random  selection  of  sightings  is  made  in  order  to  supply  conditions 
for  use  of  these  factor  tables  as  a  demonstration  of  forecasting 
and  as  a  test  of  the  reliability  of  the  data  and  of  the  forecasting 
method.  Suggestions  for  further  study  are  made. 
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1.0  Plirpoae 

The  purpose  of  this  study  Is  to  detemlne,  from  given  data, 
the  factors  affecting  the  sighting  of  surface  vessels  from  aircraft 
and  to  evaluate  the  factors  and  to  develop  a  method  of  forecasting 
the  sighting  of  surface  vessels  xinder  particular  circumstances. 
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2.0  Disouaslon  of  the  Problem 
2.1  Origin 

The  initial  impetus  toward  the  final  re8u],t8  in  this  project 
came  from  the  work  of  the  Visibility  Laboratory,  Scrlppe  Institution 
of  Oceanography,  University  of  California,  and  from  the  idea  that 
an  empirical  study  and  analysis  of  actual  sightings  of  targets  would 
be  of  assistance  in  the  Laboratory's  research  into  vision,  visibility, 
perception  and  recognition  theory.  The  first  opportunity  to  develop 
this  idea  came  with  the  discovery  (1956)  of  a  series  of  reports  of 
the  sighting  of  submarines  by  aircraft,  made  under  the  supervision 
of  Captain  Dayton  Brown,  U.S.N.R.,  Retired,  U.  S.  Navy  Electronics 
Laboratory,  A  preliminary  study  of  this  data  indicated  the 
feasibility  of  statistical  analysis,  althoxigh  there  wore  some  dis¬ 
advantages  inherent  in  the  data.  For  example  an  insufficient  number 
of  conditions  were  reported  and,  most  important,  the  data  were 
classified  under  national  security  acts.  The  classification  could 
not  be  removed.  Since  precisely  this  type  of  study  had  been  chosen 
as  a  thesis  project,  the  security  difficulty  was  frustrating.  How¬ 
ever,  enough  work  was  done  in  the  study  cited  to  demonstrate  the 
possibilities  of  extending  such  an  Investigation.  It  should  be 

ltspldcs.^i.on.5  ojT  ini'ti.H.l  yoric  ysrs  born.?  oyt 
by  the  eventual  conclusions  presented  in  this  paper. 

At  about  the  time  that  the  security  aspect  of  the  first  material 
was  becoming  rather  discouraging,  there  was  a  fortunate  development 
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(1958)  In  that  a  similar  projeot  was  reported  to  have  been  ooapletod 
recently  by  the  U.  3.  Coast  Guard  (Appendix  U  and  V).  Investigation 
showed  that  it  was  a  very  thorough  reporting  of  practically  all  of 
the  factors  affecting  visibility  of  surface  craft,  that  could  be 
readily  determined  by  a  trained  observer.  The  reports  covered 
almost  10,000  sightings  of  surface  craft  from  both  surface  and  air 
craft. 

The  results  of  this  projeot  had  been  given  to  the  Operations 
Evaluation  Group,  U.  S.  Navy,  for  analysis.  A  visit  to  the 
Group  (1958)  elicited  the  assistance  of  Dr.  J.  H.  Engel,  Deputy 
Director,  He  explained  the  scope  of  the  Group’s  work  and  vras  kind 
enough  to  turn  over  an  abstract  of  the  data  on  IBM  cards.  Frcao  his 
description  of  their  work,  it  appeared  that  more  could  be  done  than 
they  had  planned.  Thus  emerged  the  possibility  of  developing  the 
Coast  Guard  data  into  this  thesis. 

The  Operations  Evaluation  Group  has  made  an  internal  report  of 
their  study  of  the  data,  but  this  material  has  not  been  released. 

It  is  sufficient  to  say  that  the  approach  was  different  from  that 
of  the  stxidy  presented  heroin,  and  there  is  a  marked  difference 
between  results,  in  general,  and  between  interpretations  in  several 


cases. 
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2.2  Approach  to  the  Study 

The  original ^goal  of  the  study  waa  a  statistical  analysis  of 
the  various  factors  affecting  sighting  of  ships  from  aircraft  and 
deterndziation  of  the  effect  of  'tiiese  factors  on  the  sighting  range ^ 
thus  the  sighting  range  became  the  dependent  variable  in  the 
relationships  and  the  factors  affecting  it  became  the  independent 
variables.  After  consideration  of  time  and  facilities  available) 
and  the  appropriate  scope  of  the  study,  it  was  decided  to  analyse 
first  the  sighting  range  in  terms  of  each  variable  in  turn,  and 
then  to  determine  which  variables  had  significant  effects.  The 
data  would  not  be  analyzed  further  unless  required  for  a  eatlsfao- 
tory  completion  of  the  study  within  the  stated  limits.  Furthermore 
it  was  decided  that  detailed  consideration  of  tlie  Interaction  of 
subolaBBlfloatlons  of  the  variables  would  be  beyond  the  scope  of 
such  a  study.  Finally,  there  did  not  seem  to  be  a  place  for 
oonsideration  of  visibility  theory  at  this  stage  of  development. 

The  statistical  oonoopts  to  be  used  were  originally  planned 
to  be  central  tendency,  and  variance.  The  variation  of  these 
meaB^u*es  with  changes  In  each  variable  were  to  be  detennined.  The 
work  was  to  be  done  on  an  automatic  desk  calculator,  leaving  a  more 
thorough  computer  investigation  for  a  later  study  based  on  the 
results  of  this  one.  such  a  pi^limlxmry  development  had  been  made 
for  the  abandoned  Brown  data  (Section  2.1)  and  had  indicated  the 
probable  value  of  this  approach. 
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2.3  Development  of  the  Approach 

While  prelinlnary  proossslng  of  the  data  vae  proceeding  (1960)f 
there  was  a  requirenent  in  the  psyohophyeioa  research  of  the 
Vision  Branch,  Visibility  Laboratory,  for  a  computer  treatment  by 
probit  analysis  of  perceptual  threshold  experiments.  In  this 
application  a  visual  threshold  is  that  point  at  vhioh  there  is  an 
arbitrary  probability  of  seeing  a  given  target,  A  probit  analysis 
progras  (Section  2.4)  was  made  for  the  Burroughs  220  computer  based 
on  a  previously  developed  automatic  calculator  method  of  the  author. 
A  consideration  of  this  program  indicated  that  it  would  be  suitable 
for  a  study  of  the  Coast  Guard  data,  and  it  was  immediately  evident 
that  a  threshold  study  of  the  data  would  be  superior  to,  ajid  more 
useful  than,  a  straight  mean  and  standard  deviation  study.  It  would 
have  been  intellectually  uneconomical  not  to  use  the  fine  tool 
(problt  analysis)  that  was  at  hand.  So  the  standard  statistical 
approach  and  the  desk  calculator  were  relinquished  in  favor  of  pro- 
bit  analysis  on  the  high  speed  computer. 
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2.4  Probit  Analysis 

The  probit  analysis  method  comos  from  several  diverse  sources. 
Essentially  it  is  a  method  of  fitting  distribution  functions  to 
weighted  data  obtained  from  experiments.  Early  phases  of  Its 
development  were  motivated  by  the  requirement  for  the  dete-rmination 
of  kill  dosages  of  insecticides.  The  problt  method  is  based 
mathematically  on  the  maximum  likelihood  method  for  estimation  of 

« 

parameters.  The  analysis  was  developed  principally  by  Fisher  (8)» 
Garwood  (10),  and  Finney  (7),  among  others  to  be  mentioned  later. 

The  transformation  from  experimental  frequencies  to  normal  deviates 
was  introduced  by  Hazon  (11),  and  'Whipple  (24)  by  graphical  means 
and  later  developed  analytically,  as  now  used  in  problt  analysis, 
by  Wright  (25)  and  Gaddura  (9),  apparently  independently.  The 
weighting  method  involved  in  the  transformation  of  experimental 
froquenoies  is  due  to  Huller  (17)  and  was  rigorously  developed  by 
Urban  (21).  Bliss  (1  and  2)  was  responsible  for  the  name  'problt" 
and  for  a  general  description  of  the  process. 

In  the  development  of  the  probit  analysis  method,  equations 
of  estimation  wore  derived  by  the  use  of  the  maximum  likelihood 
method.  A  likelihood  function  was  eet  up  whioh  was  proportional 
to  the  product  of  the  probabilities  of  empirical  ratios  of  number  of 
responses  to  ntmber  of  presentations  of  stimuli.  Taking  the 

*  Parenthesized  numbers  following  authors*  names  refer  to  the 
corresponding  numbers  in  the  bibliography. 
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logarlthB  of  this  funotiom  sifflplifled  development  and  did  not  change 
the  oharaoter  of  the  prooeae.  The  logarithm  of  the  likelihood 
function  was  then  maxialaed  with  respect  to  the  two  unknown  para¬ 
meters  1  threshold  and  'dard  deviation!  ^  equating  tlie  partial 
derivatives  to  aero.  The  result  was  a  system  of  two  slnultameous 
equations  which  could  lead  to  determining  the  unknown  parameters 
mentioned  above.  The  aitmiltaneous  system  generally  could  not  be 
solved  by  direct  methods  but  could  be  solved  approximately  by 
iteration  after  expanding  in  Taylor's  series  and  using  trial  values 
of  the  parameters.  The  solution  of  the  eqxiatlons  was  further 
simplified  by  substituting  for  the  threshold  and  standard  deviation 
their  equivalents  in  terms  of  the  slope  and  Intercept  of  the  linear 
transformation  from  the  stimulus  domain  to  the  normal  domaini  slnoe 
the  new  trial  parameters  oould  be  found  readilyi  either  graphically! 
or  analytioally  by  the  method  of  least  squares,  or  other  fitting 
method.  The  solution  of  the  transformed  system  of  simultaneous 
equations  was  then  formulized  after  Introduotlon  of  a  device  called 
the  working  problt!  which  further  simplified  the  fomiulatlon. 

The  probit  analysis  method  is  adaptable  to  many  experimental 
problems  involving  cumulative  quantal  data.  The  method  Is  in  use 
In  payohophysioal  research  at  the  University  of  Michigan  and 
the  Unlvaralty  of  Ralifomia.  It  has  been  adapted  for  desk  calcu¬ 
lator  use  by  Kincaid  and  Blackwell  (13)!  '^y  the  author  (18)} 

to  digital  computers  by  Mol.dauer  and  Klnoald  (16)  for  MIDAC!  by 
Lamphlear  and  Wendel  (14)  for  IBM  650!  and  by  the  author  for 
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BurroughB  220  (future  publioation)  and  CDC  1604  ai\d  IBM  7090 
(future  publioation). 

The  adaptation  used  in  this  project  is  that  of  the  author  for 
the  Burroughs  220  (future  publication)  which  in  turn  is  based  on 
his  adaptation  for  desk  oaloulatoro  (18).  Olven  cumulative  poeitive 
responses  to  stimuli  over  the  range  of  the  stimuli,  ratios  or 
empirical  probabllltleB  are  calculated  and  transformed  to  abscissas 
of  the  normal  distribution  function,  normal  (0,1),  that  is  with 
mean  of  zero  and  standard  deviation  of  one.  Finney  (7)  uses  a 
5-biased  normal  distribution,  normal  (5,1),  to  avoid  negative 
abscissas.  The  resulting  abscissa  is  called  the  probit,  hence  the 
name  of  the  method.  In  the  digital  computer  treatment  negative 
abscissas  are  no  disadvantage  so  the  5-bia8  is  not  used,  though  the 
liberty  is  taken  of  calling  this  abscissa  a  probit  also,  since  the 
use  is  the  same.  The  probit  which  has  been  determined  is  termed 
the  empirical  probit,  A  trial  linear  tranaformatlon  function  of 
the  form  y  =  a  +  bx  to  transform  from  the  experimental  domain  to 
the  normal  domain  is  determined,  in  this  mettiod,  by  the  method  of 
least  squares.  The  resulting  probits,  y  ,  corresponding  to  the 
experimental  stimulus  points,  x  ,  are  termed  trial  problts.  Having 
determined  the  trial  probits  for  the  range  of  stiinull,  the  weighting 
factor  is  applied  to  account  empirical 

probabilities  toward  tho  tails  of  the  distribution  function.  All  of 
the  requirements  are  now  available  for  the  solution  of  the  maximum 
likelihood  system.  The  solution  is  facilitated  by  the  introduction 
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of  tho  working  probity  sonetinaB  and  inoorreotly  called  the  corrected 

problti  which  allows  a  simply  calculated,  next  approximation  to  the 

parameters  of  the  transformation  function.  Approximate  values  of 

the  threshold  and  standard  deviation  are  determined  directly  from 

2 

the  parameters  and  are  used  to  make  a  y  relation 

between  the  empirioal  frequency  ratio  and  the  ratio  calculated  from 
the  threshold  and  standard  deviation.  If  the  teat  shows  a  signifi¬ 
cant  difference,  the  transformation  parameters  are  used  as  corrected 
trial  parameters  to  reenter  the  process  in  order  to  improve  the 

results  by  iteration.  The  computer  application  used  in  this  project 

2 

oontinuee  iteration  until  an  acceptable  y  value  results  or  until 

2 

the  process  begins  to  diverge.  If  an  acceptable  y  value  results, 
the  final  approximation  of  the  thi'eshold  and  standard  deviation  are 
recorded. 

The  flow  of  the  computation  is  shown  in  Appendix  W. 

It  may  be  of  interest  that  the  computer  process  requires  a 
computation  time  of  2  to  3  seconds,  since  card  input  of  data  is 
folded  into  the  oomputatlon.  Print-out  time  is  about  30  seconds 
since  a  Soroban  teletype  printer  is  used.  Use  of  an  IBM  407 
printer  would  cut  overall  problem  time  to  about  5  seconds. 
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3*0  Dlsouaaion  of  the  Analyais  Prooeaa 
3.1  Sorting 

The  initial  step  in  processing  the  mass  of  data  was  to  separata 
the  air  sightings  from  the  10,000  cards  that  included  both  air  and 
surface  sightings.  This  was  done  on  an  IBM  mechanical  card  sorter 
and  the  cards  were  sorted  both  by  type  of  observing  vehicle  and  by 
observer  to  assure  that  no  surface  sightings  were  included  in  the 
3,465  cards,  each  of  which  documented  an  air  sighting.  This  sorting 
was  done  before  the  decision  was  made  to  use  probit  analysis.  The 
data  were  next  sorted  to  classify  them  for  a  standard  statistical 
analysis.  This  sort  was  with  respect  to  the  dependent  variable, 
sighting  range.  The  range  of  the  variable  was  checked  and  it 
appeared  that  the  last  class  should  include  22  miles  and  more,  since 
any  further  classes  would  contain  too  few  sightings  for  reliable 
analyais.  This  procedure  was  very  convenient:  the  optimum  number 
of  classes  from  the  standpoint  of  both  reliability  and  practicality 
is  usually  taken  as  from  ten  to  fifteen;  thus,  using  an  increment  of 
two  miles  from  zero  to  twenty-two  miles,  there  were  twelve  classes. 
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3t2  Counting  and  Tabulating 

The  oard  sorter  was  equipped  with  a  counter  and  the  next  proce¬ 
dure  was  to  sort  each  class  of  the  dependent  variable  Into  olassas  of 
an  Independent  variable  and  then  rerun  each  class  of  the  Independent 
variable  to  determine  the  count.  Here  the  number  of  olasses  does 
not  depend  on  statistical  theory  but  on  a  logical  division  of  the 
ranges  of  the  Independent  variable  so  as  to  give  useful  results  in 
the  case  of  a  continuous  type  variable.  Some  of  the  Independent 
variables  were  inherently  broken  down  into  discrete  classes  as  in 
the  ease  of  the  visual  aids  used.  The  result  of  this  sort  and  count 
on  the  seventeen  independent  variables  was  tabulated  and  produced 
seventeen  frequency  matrices. 
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3i3  Arranging  for  Problt  Analyais 

One  more  process ing  step  remained  before  the  data  could  be  pre¬ 
sented  to  the  computer.  Problt  analysis  treats  cumulative  distribu¬ 
tion  functions  and  not  Gaussian  frequency  functions.  A  distribution 
matrix  was  oonpllcd  with  each  class  of  the  independent  variable 
becoming  a  distribution  function  vector  with  respect  to  the  sighting 
range.  In  other  words,  the  sighting  range  then  became  the  indepen¬ 
dent  variable  for  problt  pui^ioses,  and  the  frequency  within  the  class 
of  the  independent  variable  became  the  dependent  variable,  or  distri¬ 
bution  function  of  the  sighting  range.  This  resulted  in  147  probit 
problems. 

In  sorting  and  counting  on  the  independent  variables,  cards 
coded  for  no-report-entry  or  anomalous  entry  were  tabulated  and 
given  a  balancing  check  against  the  over-all  count.  They  were  also 
included  in  the  probit  analysis  of  the  variable  as  an  extra  class, 
but  no  significance  was  noted  other  than  that,  as  sighting  conditions 
become  more  difficult,  obsejrvers  tend  to  be  more  meticulous  (which 
was  Inferred  from  the  lower  proportion  of  faulty  cards) . 

The  data  were  ready  for  computer  manipulation  at  this  point. 
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3.4  Computer  Prooesslng  and  Probleras  in  the  Problt  Analysis 

The  data  ware  now  punched  into  cards  for  entry  to  the  computer. 
The  original  olasslfleation  was  followed,  each  sighting  range  incre¬ 
ment  of  two  miles  becoming  a  point  determining  the  distribution.  In 
general  this  gave  twelve  points  on  the  abscissa  for  even  numbered 
miles.  The  distribution  functions  were  found  to  be  of  the  log-normal 
type,  considered  to  be  a  result  of  the  exponential  attenuation  of 
radiance  and  contrast.  The  computer  program  provided  for  both  nor¬ 
mal  and  log-normal  distributions.  The  input  data  called  for  a 
selection  of  log -normal  analysis  and  the  computation  was  started 
on  this  basis.  Apparently  a  major  task  in  the  processing  had  been 
completed  and  the  results  of  the  computation  were  ready  for  analysis 
of  effects  and  trends.  Unfortunately  this  proved  not  to  be  the  ease. 

The  computer  output  consists  of  the  parameters  of  the  problt 

transfer  function,  the  threshold,  the  fiducial  limits,  the  standard 

deviations  of  the  data,  threshold,  standard  deviation  and  parameters, 

2 

and  last,  but  far  from  least,  the  x  measure  of  goodness  of  fit. 

2 

An  examination  of  the  x  measure  showed  that  only  a  little 
over  half  of  the  functions  were  fitted  at  the  0,05  level,  which  had 
been  selected  as  the  acceptable  criterion.  The  unacceptable  func¬ 
tions  were  cheeked  and  end  points  in  the  tails  of  the  distributions, 
that  showed  very  small  numbers  in  the  frequency  table,  were  stripped 
out.  This  is  usual  procedure  In  probit  analysis,  since  it  does  not 
affect  the  res\ilts,  and  the  instability  in  the  tails  may  affect  the 
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X  measure,  not  the  essential  character  of  the  distribution. 

Reruns  were  made  of  the  unacoeptable  data  sets  and,  while  there  was 
a  marked  gain  in  production,  there  was  still  about  a  quarter  of  the 
sets  that  were  not  acceptable. 

It  was  at  this  point  that  what  had  been  noticed  as  an  interest¬ 
ing  sidelight  in  the  data  counting  became  a  matter  of  crucial  impor¬ 
tance  j  it  had  been  evident  in  the  sorting  and  counting  that 
observers  tend  to  estimate,  or  rotind  off,  to  multiples  of  five,  and 
the  stacks  in  these  bins  were  disproportionately  high.  It  was  now 
evident  that,  with  the  relatively  fine  definition  of  a  two  mile 
increment  in  sighting  range,  this  tendency  was  introducing  an  extra¬ 
neous  scallop  in  the  distribution  functicas. 

An  obvious  method  of  removing  this  soallop  was  to  increase  the 
increment  in  the  sighting  range  so  that  the  data  divisions  vere 
located  at  points  where  the  data  character  was  consistent,  such  as 
on  multiples  of  five  miles.  Here  the  sorting  on  even  miles  might 
have  been  a  possible  handicap,  for  the  only  alternative  to  resorting 
and  recounting  the  basic  carxl  deck  was  to  choose  the  divisions  of 
the  existing  sort  that  included  the  multiples  of  five  miles.  This 
course  was  decided  for  trial,  and  all  of  the  data  sets  were  set  up 
again  with  abscissa  points  at  zero,  four,  ten,  fourteen  and  twenty 


miles,  since  these  ixioludeu  Uie  five  mile  divisions.  The  results  of 


this  run  woro  gratifying,  for  only  about  a  tenth  of  the  problems  were 


not  acceptable  at  the  0.05  level  and  most  of  these  were  acceptable  at 


the  0.01  level.  This  last  level  was  not  considered  acceptable  and 
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the  tall  points  with  very  few  numbers  were  stripped  out.  Most  of 
these  were  equivalent  to  probabilities  of  less  than  0.01  and  in  no 
case  as  muoh  as  0.03.  It  la  of  Interest  to  note  here  that  taany  of 
the  acceptable  problems  Included  probabilities  of  less  than  0.01, 
which  is  not  usually  expected.  A  rerun  was  made  of  the  stripped 
problems  and,  out  of  those  that  might  be  expected  to  give  good 
results,  only  three  remained  unacceptable  at  the  0.05  level.  That 
la,  the  few  others  remaining  had  less  than  six  sightings  to  a  pro¬ 
blem. 

While  one  can  solve  for  any  probability  threshold  in  probit 
analysis,  depending  on  the  needs  of  the  analysis,  in  this  study  the 
0.5  threshold  is  found.  This  is  also  known  as  the  "505?  threshold," 
or  the  "mean  threshold,"  It  is  the  point  at  which  the  probability 
of  sighting  is  0.5,  or  that  point  at  which  an  observer  is  as  likely 
to  see  as  not  to  see  an  object.  The  conversion  of  the  0.5  threshold 
to  a  threshold  of  any  other  probability  is  given  later. 
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3.5  DevBlopment  of  the  Produced  Data 
3.5.1  Tables 

The  computer-proceaaed  results  are  tabulated  ( Appendices 
A  through  P)  to  ahow  the  Initial  data  and  the  probit  analysis  para¬ 
meters:  threshold  (T),  standard  deviation  of  the  distribution  (s), 
standard  deviation  of  the  threshold  (s,p)  and  the  chi-square  measure 
of  goodness  of  fit  (x  )•  Where  a  small  probability  tail  value  is 
dropped,  it  is  shown  in  parenthesis. 

The  mean  standard  deviation  (S) ,  the  standard  deviation  of  the 

standard  deviations  (a  ),  the  standard  deviation  of  the  mean  stan- 

s 

dard  deviation  (s-)  were  also  added  to  the  tables. 

'  3 
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3.5 <2  Graphical  Display 

The  thresholds  are  plotted  on  graph  paperi  in  a  linear  plot 
or  in  a  logarithmic  plot  if  this  appears  appropriate.  (Appendices 
B  through  P) .  This  procedure  aids  in  determining  the  type  of  empiri¬ 
cal  function  to  be  fitted  to  the  data.  Where  a  sine  or  cosine  type 
curve  appears  suitable,  double  plotting  is  used  to  facilitate  visual 
interpretation.  In  this  ease  the  original  plot  is  made  and  the 
supplement  or  complement  of  the  angle  is  used  to  plot  the  correspond¬ 
ing  thresholds  on  the  same  abscissa  axis. 
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3 •5.3  Fitting  Empirleel  Funotions 

After  a  survey  of  the  data  graphs,  empirical  funotions  are 
fitted  to  the  threshold  data  using  standard  least  square  methods  of 
fit.  In  the  case  of  parabolio  functions,  Cholesky's  method  (19)  of 
reduction  of  mntrioea  is  used  in  the  solution  of  the  normal  least 
square  equations.  This  method  greatly  facilitates  the  solution  of 
simultaneous  linear  systems.  It  is  described  by  Salvador!  and  Baron 
(19).  The  simplest  function  is  chosen  which  gives  a  high  correlation 
coefficient  and  preserves  the  character  of  the  data.  These  functions 
are  added  to  the  tables  along  with  a  factor  function  (f)  which  is  the 
empirical  function  normaliaed  to  the  normal  range  which  will  be  des¬ 
cribed  in  the  section  on  the  general  distribution  of  sightings 
(Section  4.1).  Functional  threshold  values  are  then  calculated  and 
added  to  the  graphs  for  comparison  purposes. 
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3.5.4  Testing  for  Slgnlficanoe 

Maximum  and  minimum  functional  threshold  values  are  tested  to 
determine  whether  or  not  the  effect  of  the  variable  is  significant. 
The  t  test  for  difference  of  two  means  is  used,  based  on  Student's 
t  distribution.  The  only  case  of  insignificant  effect  is  that  of 
wind  azimuth  which  is  to  be  discussed  In  detail  in  that  section. 
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3.5 *5  Notes  on  Results 

It  is  very  Important  to  keep  In  mind  that,  while  the  threshold 

T  is  in  the  linear  domain  and  dimensioned  in  nautical  miles, 

3  ,  Bj  ,  a  ,  Sg  are  in  the  logarithmic  domain  and  are  dlmenslonleBS 

2 

constants  in  this  application.  ^  ,  t  and  F  also  apply  to  the 

logarithmic  domain. 

To  use  the  standard  deviation  as  a  probability  tool  In  deter¬ 
mining  a  threshold  of  other  than  0.5,  the  following  relationship 
holds: 

T(p)  =  T, 

where  T(p)  is  the  threshold  for  the  desired  cumulative 
probability  p  . 

T  is  the  threshold  which  is  always  for  a 
p  =  0.5  in  this  study. 

D(p)  is  the  normal  deviate  for  the  probability  p  . 
s  is  the  standard  deviation. 

The  reason  for  the  negative  exponent  in  the  above  relationship  la 
that  these  are  reverse  distribution  functions;  the  distribution 
function  is  monotonically  decreasing  in  the  positive  abscissa 
direction. 

The  above  formula  is  specialized  and,  if  the  user  is  working  in 
other  units  than  nautical  miles,  T  must  be  converted  to  the  new 
units  in  order  to  compute  T(p).  The  e  la  not  converted  in  any  way. 
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It  1b  interesting  to  note  the  extreme  sensitivity  of  the  problt 

2 

analysis  technique  in  conjunction  with  the  x  test.  Twenty  of  the 
problemsi  not  acceptable  on  the  first  run,  are  made  acceptable  by 
removing  the  very  small  numbers  in  the  extreme  tails  of  the  distribu¬ 
tions.  In  nineteen  of  the  problems,  this  change  results  in  an 

2 

average  reduction  of  x  from  10.253  to  1.894  while  there  is  an 

average  absolute  ohange  of  only  0,00965  of  the  value  of  the  threshold. 

This  result  definitely  indicates  the  inherent  stability  and  excellence 

of  the  method.  In  the  one  remaining  problem,  removing  the  small  tail 

2 

number  proportionally  changes  the  threshold  0.186  while  reducing  x 
from  72,157  to  the  magnitude  of  10~^^.  It  is  surely  a  coincidence 
with  the  criterion  of  acceptance  that  this  one  reading  is  0.05  of  the 
twenty  readings  In  this  category. 
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4.0  The  Analysis;  Discussion  of  Graphs  and  Tables 

In  this  section  a  detailed  discussion  is  made,  where  appropriate, 
of  the  input  data  to,  and  results  of,  the  probit  analysis.  In  par¬ 
ticular,  comments  are  made  on  the  graphical  evidence  and  empirical 
function  determination.  The  empirical  threshold  function,  T(.);  the 
normalized  threshold  or  factor  function,  f(.);  and  the  correlation 
coefficient,  r  are  stated  in  this  section.  Any  other  significant 
observations  are  included  with  each  variable  treated. 

The  mean  of  the  computed  standard  deviations  is  entered  in  the 
appendix  tables  for  information  and  possible  use,  since  the  standard 
deviations  of  the  thresholds  appear  generally  level  enough  so  that 
the  mean  standard  deviation  might  well  be  used  through  the  range  of 


most  of  the  variables, 
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4>1  General  Distribution  (See  Appendix  A) 

In  this  case  all  reoorde  that  have  sighting  ranges  are  lumped 
together  and  a  distribution  set  up.  The  data  shows  a  typically  log¬ 
normal  oharaoter.  The  problt  calculation  gives  a  threshold  of  6.599 

2 

nautical  miles  with  a  x  4.308,  which  is  acceptable  at  the  0.05 
level. 

This  threshold  for  all  sightings  is  based  on  3»465  sightings 
taken  on  the  Atlantic  and  Pacific  coasts  from  Puerto  Rico  to  Alaska 
and  should  be  a  representative  measure  of  central  tendency  of  the 
class  of  all  sightings  of  surface  vessels  from  aircraft.  This  con¬ 
clusion  is  borne  out  by  the  estimate  of  the  coefficient  of  variation 
of  the  threshold,  0.0015.  Hence  this  threshold  range,  6,599  miles. 
Is  taken  as  the  standard  reference  range  in  this  study  and  is  used 
as  the  normalising  range  where  normalization  is  carried  out,  as  in 
the  development  of  the  threshold  factors.  This  range  is  called  the 
"normal  range,"  ,  end  is  so  marked  on  the  graphs. 
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4.2  Meteorologioal  Vlelbility  (Appendix  B) 

Meteorological  vlalblllty  Is  defined  as  follows  by  the  U.  S, 
Weather  Bureau  (23) i  "...  the  greatest  distance  toward  the  horizon 
that  prominent  objeots  euoh  as  mountains,  buildings,  towers,  eto., 
can  be  seen  and  identified  by  the  normal  eye,  unaided  by  special 
optical  devices,  such  as  binooulare,  telescopes,  glare  ellrainators, 
goggles,  etc.,  and  which  distance  must  prevail  over  the  range  of  more 
than  half  the  horizon."  Granting  that  this  procedure  does  not  present 
a  rigorous  measure  of  atmospheric  turbidity,  it  is  still  the  best 
one  now  in  general  use  and  is  the  one  in  which  all  observers  are 
trained,  have  used,  and  that  most  will  use  for  some  time  to  come. 
However  it  should  be  noted  that  precise  definitions  of  visibility  have 
been  formulated  in  terms  of  "meteorological  range,"  a  quantity  that 
can  be  defined  analytically,  measured  instrumentally  and  used  in  the 
solution  of  visibility  problems. 

The  0-1  mile  visibility  classification  has  only  three  sight¬ 
ings,  one  at  six  miles  sighting  range  and  two  at  zero  miles.  While 
the  cumulative  data  are  routinely  introduced  in  the  probit  analysis, 

this  classification,  as  might  be  expected,  gives  a  T  of  the  order 
—11 

of  10  ,  s  of  10  ,  8^  of  10  .  Any  further  consideration  of 

this  classification  is  ruled  out.  There  is  an  apparent  anomaly  in 
the  data  for  the  classification  of  10  -  19  miles  vlBiblllty  between 
the  sighting  ranges  of  ten  and  twenty  miles.  Numerous  efforts  with 
various  combinations  of  input  result  in  no  convergence  at  the  0.05 
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level.  The  count  has  been  rechecked  and  the  only  possible  conclu¬ 
sion  is  that  this  ia  one  of  the  thing’s  that  sometimes  happens.  This 
la  the  only  case  in  this  study  where  the  0.05  aoaeptanoe  criterion  is 

relaxed.  The  data  inaluded  for  this  olaesifioation  in  the  appendix 

2 

are  acceptable  at  the  0.01  level  with  ax  10.094  and  are  the 

2 

best  fit  obtained.  The  next  best  fit  ia  at  a  x  15.337  which  has 
a  proportional  difference  of  0.00239  from  the  accepted  threshold. 
Since  the  stability  of  the  mean  is  established  thus,  and  because  of 
the  fact  that  all  other  results  group  very  closely  about  the  accepted 
threshold,  the  number  that  appears  in  the  table  is  accepted.  This 
procedure  is  admittedly  subjective  reasoning  to  a  certain  extent  but 
it  ig,  a  considered  judgment. 

A  plot  on  semi-logarithmic  paper  shows  a  distinct  linearity  of 
the  type  T(V)  =  a  +  b  InV.  A  least  square  fit  gives  this 

equation  I 

T(V)  =  3.476  InV  -  2.064,  V  =  Visibility 

f(V)  =  0.527  InV  -  0.313 

r  =  0.969  ,  in  the  logarithmic  domain. 

The  sighting  range  may  well  approach  an  asymptote  as  visibility 
increases,  but  there  is  insufficient  evidence  to  determine  this 
effect  at  this  time. 

A  point  of  Interest  in  the  threshold  table  is  the  sharp  and 
(at  first  thought)  surprising  drop  in  sighting  range  with  unlimited 
visibility.  It  would  appear  that  when  an  obseinrer  has  no  definite 
landmarks  but  good,  clear  air,  he  terms  this  condition  "unlimited," 


28 


SIO  Reft  62-13 


Howevaff  If  ha  states  visibility  is  70  milesi  he  doubtlesa  identifies 
an  objeot  at  that  knovm  range.  May  we  assuine  that  the  olaasifioation 
'*\inlimited'’  really  means,  "I  think  I  can  see  a  long  way"?  Meteoro- 
logikts  agree  with  this  assumption  and  consider  that  a  report  of 
unlimited  visibility  should  be  taken  to  mean  "more  than  fifteen 


miles." 
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4.3  Altltvide  (See  Appendix  C) 

The  data  on  thle  variable  require  diecuseion  only  in  the  9»000 
foot  claBsification.  No  difficulty  is  involved  in  the  eolution  of 
any  of  the  other  olaseea.  The  usual  run  at  9»000  feet  including  the 
nultiples  of  five  miles  sighting  range  gives  anomalous  results.  Pre¬ 
viously!  increments  of  two  miles !  aooeptable  results  had 

2 

appeared,  and  even  better  results  with  respect  to  x  ooourred 

on  increments  of  four  miles.  This  last  interval  is  the  number  used. 
The  irregularity  results  no  doubt  from  only  seven  sightings  at  this 
altitude.  The  data  shows  a  linear  trend  In  a  seni-logarlthmlo  domain 
with  a  funeticnal  form  of  T(A)  =  ae'’^.  The  best  fit  function  Isi 

T(A)  =  5.3B5  ,  A  =  altitude  in  lOOO's  of  feet. 

f(A)  =  0.816 

r  =  0.963,  in  the  fitting  domain. 

There  is  no  apparent  reason  for  the  reverse  trend  from  2,000 
to  4,000  feet,  such  as  a  general  hemispheric  tendency  to  a  3,000  foot 
haze  layer.  The  linear  plot  shows  close  grouping  about  the  fitted 
line.  Hence  no  attempt  is  made  to  fit  a  cubic  equation.  This  might 
be  a  natter  for  further  Investigation. 
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4«4  Ship  Siae  (Sae  Appendix  D) 

The  data  here  appear  unremarkable,  and  the  probit  analysis  pre¬ 
sents  no  difficulty,  requiring  reruns  only  where  tliere  are  four  and 
five  tail  sightings  in  the  classes:  less  than  30  feet  bright- 
colored;  30  to  60  feat  bright  colored. 

It  la  in  this  variable  that  the  only  major  criticism  of  the 
report  form  (Appendix  V)  arises.  The  dimension  of  the  variable 
changes  at  100  feet  from  length  to  tonnage.  This  discontinuity 
requires  a  transformation  from  tonnage  to  length.  With  the  help 
of  Fahey's  catalog  of  U.  S.  Navy  auxiliaries  (6),  a  fxmotlon  la 
developed  to  effect  this  transformation.  It  appears  that  the  list  of 
axuclliarles  is  a  representative  cross-section  of  ships  that  a  Coast 
Guard  patrol  would  sight.  The  function  adopted  is: 

L  =  21.4  T,y^  -  16 
w 

where  Ty  is  full  load  tonnage  and  L  is  length  in  feet.  This 
function  has  a  correlation  coefficient  of  0.985.  (See  Appendix 
D-3  and  4) 

After  transforming  the  data  tonnages  to  length,  the  threshold 
data  show  a  very  decided  linearity  in  the  semi-logarlthmio  domain 
of  the  type  T(L)  =  a  +  b  InL.  The  best  fit  function  is: 

T{L)  1.844  InL  -  1.100,  L  ■=  Length  in  feat 


f(L)  =  0.280  InL  -  0.167 


I 


SIO  Refi  6a-13  31 

r  a  0.982 

The  effect  of  shading  of  target  might  be  handled  In  a  number  of 
ways,  but  various  trials  show  that  a  simple  and  satisfactory  correc¬ 
tion  is  to  add  0.279  miles  to  the  threshold  for  bright  vessels  and 
subtract  0.279  miles  for  dark  vessels,  both  types  under  100  feet. 

There  is  insufficient  basis  to  apply  this  correction  to  larger 
targets.  The  correction  is  based  on  a  fit  made  to  the  separate 
bright  and  dark  classes. 
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4.5  Height  of  Major  Swelle  (See  Appendix  E) 

There  le  no  need  to  comment  on  the  data  here  other  than  to  men¬ 
tion  that  mean  swell  height  for  the  ten-feet-and-more  classification 
is  15.381  feet  and  is  so  used  in  the  fitting  process.  The  probit 
analysis  turns  out  well. 

The  data  show  a  euffieient  linear  trend  to  adopt  an  empirical 
function  of  the  T(S)  =  a  +  bS  type.  The  fit  produced  Isj 

T(S)  =  6.170  +  0.239  S  ,  S  =  Swell  height  in  feet 

f(S)  =  0.935  +  0.0362  S 

r  =  0,707 

The  correlation  coefficient  is  not  as  large  as  would  be  desired^ 
due  mainly  to  the  two  highest  values.  Including  those  values  might 
be  questioned,  in  view  of  the  few  sightings  in  these  olasslficationa. 
However  they  show  no  anomalies,  have  p(x  )  values  of  greater  than 
0.9  and  0.4  respectively,  and  they  straddle  the  fitted  line.  It  la 
considered  better  to  retain  them,  lacking  further  evidence. 

It  may  seem  surprising  that  one  sees  objects  better  in  higher 
swells,  but  there  may  be  three  causes  for  this.  One  is  the  relative 
motion  of  the  object  with  respect  to  the  water  masses,  which  attracts 
attention.  Another  is  that  the  higher  swells  break  up  the  grazing 
reflection  of  the  brighter  horizon  sky.  A  third  perhaps  is  that  the 
object  is  seen  against  a  surface  sloped  towaird  the  observer  which  may 
have  some  effect  by  giving  a  generally  darker  background  rather  than 
the  horizon  sky  reflection.  This  same  trend  has  been  noticed  by 


Coast  Guard  observers 
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4.6  Cloud  Cover  (See  Appendix  F) 

No  dlffloultF  ia  enoountered  with  thla  data  or  with  the  problt 

analysis.  An  inspection  of  the  linear  plot  clearly  indicates  a 

2 

parabolic  fit  of  the  form  T(C)  =  a  +  bC  +  cC  .  The  least  square 


method  gives: 

T(C) 

=  7.069  +  2.371C  -  4.717C^,  C  =  Cloud  cover 

in  decimal  fraction 

f(c) 

=  1.071  +  0.435C  -  0.715C^ 

r 

=  0.903,  in  the  linear  domain 

The  maximum  value  of  this  function,  T(0,304)  =  7.506  miles 

threshold,  bears  out  informal  Coast  Guard  impressions  that  an  obser¬ 
ver  sees  best  with  about  one-third  oloud  cover.  This  may  be  a  result 
of  diminution  of  surface  glare  and  glitter  with  enough  direct  light¬ 
ing  remaining  to  give  good  contrasts. 
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4.7  Wind  Velocity  fSee  Appendix  G) 

The  one  sighting  of  45  -  49  knots  Is  disregarded.  At  50  knots 
and  greater  there  are  only  four  sightings.  The  probit  analysis  is 
good  until  the  50  knot  results  are  inspected  carefully.  Here  the 
parameters  of  the  probit  transfer  equation  are  of  an  entirely  differ¬ 
ent  magnitude  and  character  from  the  rest  of  the  family.  Since  there 
are  so  few  sightings  (and  this  is  an  extremely  erratic  threshold  of 
low  dependabilltyj  =  0.43  in  the  log  domain)  it  is  rejected. 

The  other  two  erratic  points  are  accepted.  There  Is  no  clear  reason 
for  their  rejection. 

The  accepted  thresholds,  whan  plotted,  Indicate  a  parabolic 
trend  of  the  type  T(WV)  =  a  +  bWV  +  o(WV)^  .  The  beat  fit  is: 

T(WV)  =  '}.W  -  0.142  W  -  0.002673(WV)^ 

WV  =  wind  velocity  in  knots 

f(WV)  =  0.832  +  0,0214  WV  -  0.000405(WV)^ 

r  =  0.531 

Here  the  correlation  coefficient  is  lower  than  desirable,  due  to  the 
two  erratic  points  included.  However,  the  empirical  function  is  the 
only  simple  function  typical  of  the  trend  of  the  data.  The  maximum 
value  of  the  function,  T(26.6)  =  7.378  miles  indicates  the  best 

seeing  is  where  the  meteorologists  say  a  breeze  becomes  a  gale 
between  27  and  28  knots,  or  between  force  six  and  force  seven  winds 
on  the  Beaufort  scale.  This  conclusion  is  reasonable  since  in  this 
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veloeity  range  the  surface  beeotnes  very  disturbed  and  spray  and 
foam  appear.  This  condition  explains  the  different  character  of 
the  variable,  WV  ,  from  that  of  swell  height. 
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4.8  Relative  Bearing  of  the  Target  (See  Appendix  H) 

This  variable  is  also  termed  "eloek  code"  in  the  report  form 
(See  Appendix  V).  The  data  show  very  few  sightings  in  the  I’ear  field 
of  view.  This  sighting  lack  raises  the  question  of  the  practicality 
of  designing  patrol  planett  with  tail  observations  a  major  factor. 

The  problt  analysis  gives  good  results  although  the  character  of 
the  few  sightings  at  seven  o'clock  defeats  efforts  to  determine  an 
acceptable  threshold.  The  low  number  of  sightings  to  the  rear  and 
their  somewhat  erratic  behavior  make  it  tempting  to  discard  them. 
However,  double  plotting  shows  tee  data  are  consistent  in  character 
and  follow  a  cyclic  pattern  of  tee  type  T(B)  =  a  +  b  cos  B. 

The  best  fit  function  is: 

T(B)  =  6.046  +  0.808  cos  B 

B  =  Relative  bearing  or  30 ’(clock  code)  degrees 

f(B)  =  0.916  +  0.122  COB  B 

r  =  0.697 

Again,  the  relatively  low  r  la  the  result  of  the  erratic  rear 
sightings.  However,  a  polar  plot  of  this  function  shows  an  appa¬ 
rently  less  erratic  trend  to  the  rear. 
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4.9  Sun  Altitude  (See  Appendix  I) 

Here  the  only  case  demanding  connnent  is  that  at  90  degrees  where 
there  are  only  six  sightings.  The  probit  analysis  does  not  give  an 
acceptable  fit  and  this  classification  is  disregarded.  Otherwise  the 
probit  analysis  is  normal. 

A  survey  of  the  plotted  thresholds  indicates  an  empirical  func¬ 
tion  of  the  type  T(SA)  =  a  +  b  sin  (SA),  and  the  least  square 
method  produces; 

T(SA)  =  7.795  -  1.564  sin  (SA),  SA  =  Sun  altitude: 

0°  Z  SA  Z  90° 

f(SA)  =  1.208  -  0.238  sin  (SA) 

r  =  0.915  in  the  fitting  domain 

Thus,  it  appears  that  seeing  deteriorates  with  increasing  sun  alti¬ 
tude.  This  conclusion  seems  reasonable  since  a  low  sun  results  in 
stronger  Internal  contrasts  in  the  target  while  a  high  sun  gives  a 
flatter  lighting. 
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4.10  Relative  Bearing  of  Sun  (See  Appendix  J) 

This  bearing  is  with  respect  to  the  target.  The  data  hare  is 
worthy  of  comment  in  that  the  total  number  of  sightings  in  each 
olassification  of  bearing  alternate  in  magnitude  with  the  high  number 
applying  to  the  class if icatione  containing  the  four  major  divisions  of 
the  circle,  multiples  of  90  degrees.  There  is  no  probit  trouble  and 
every  problem  develops  acceptably  the  first  time. 

Double  plotting  of  the  data  (as  in  Section  3.5,2)  clearly  indi- 
oates  an  empirical  oyelic  function  of  the  type  T(SB)  =»  a  + 

+  oo8(iT  -  SB).  The  oharaoter  of  the  data  necessitates  the  phase 
shift.  The  fitting  process  gives: 

T(SB)  =  7.043  +  0.342  oo3(n  -  SB),  SB  =  Sunbearing  in 

degrees 

f(SB)  =  1.067  +  0.0518(n  -  SB) 

r  =  0.987  in  the  fitting  domain 

This  conclusion  shows,  as  might  bo  expected,  that  an  observer  sees 
better  with  the  sun  behind  him  as  he  looks  at  the  target.  However, 
the  change  from  minimum  to  maximum  of  6.701  miles  to  7.385  is  not  as 
great  as  one  might  expect. 

The  question  may  be  raised  why  the  minimum  threshold  la  greater 
than  the  normal  range.  Does  one  expect,  on  the  basis  of  this  beha¬ 
vior,  always  to  see  better  than  the  normal  range?  Examination  of  the 
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data  leads  to  the  explanation  that  in  total  overoaet,  oloud  oover 
of  1.0,  relative  bearing  of  eun  with  reepeot  to  target  la  frequently 
not  reported  and  under  these  overoaet  conditions  the  threshold  Is 
much  less  than  the  normal  range. 

A  further  point  of  Interest  Is  that  the  alternate  claeslfloations 
with  low  numbers  of  sightings,  mentioned  above  in  the  first  para¬ 
graph,  are  also  the  class If loatlons  having  the  larger  thresholds. 
Perhaps  this  distribution  means  that  those  who  are  meticulous  In 
noting  exact  bearing,  rather  than  the  nearest  90  degree  direction, 
are  also  more  meticulous  and  alert  in  their  search  operation. 
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4.11  Proportional  Waka  Size  (See  Appendix  K) 

This  le  the  last  of  the  conditions  that  lead  to  a  continuous 
type  threshold  function  and,  in  effect,  the  least  important.  The 
choice  of  proportional  instead  of  actual  wake  size  was  due,  doubtless, 
to  plan  rather  than  fortune.  This  proportional  iseasure  automatically 
rules  out  such  aberrations  as  might  result  from  large  slow  ships 
with  little  wake  in  comparison  with  small  fast  vessels  leaving  large 
wakes.  The  effect  here  la  small  though  significant  statistically. 

The  data  appears  good  as  is  proven  by  no  failures  on  the  first 
problt  runs.  A  survey  of  the  plotted  thresholds  indicates  an  expo¬ 
nential  curve  approaching  a  non-zero  asymptote.  The  type  cuirve  is 
T(WS)  =  a  +  .  Taking  the  asymptote,  T(WS)  =  a  ,  to  be  the 

threshold  of  the  classification  defined  as  greater  than  twice  the 
length  of  the  vessel,  the  following  function  develops i 

T(WS)  =  7.295  -  1.066e"^’^^^^'^®^  WS  =  Wake  size, 
as  a  multiple  of  ship  length 

f(WS)  =  1.105  - 

r  =  0.920  in  the  fitting  domain 

The  correlation  Is  even  better  in  the  linear  domain  so  no  attempt 
is  made  to  determine  the  asymptote  analytically,  a  doubtful  procedure 
at  best  with  as  few  data  points  as  are  available  here. 
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4.12  Wind  Azimuth  (See  Appendix  L) 

It  Is  unfortunate  that  the  definition  of  wind  direction  as  an 
azimuth  was  prescribed  In  the  report  form  (see  Appendix  V) .  It 
seems  very  unlikely  that  wind  azimuth  would  ehow  a  significant  effect 
on  threeholde  of  sightings,  whose  azimuths  must  be  assumed  to  be  more 
or  less  random.  A  much  batter  and  more  valuable  definition  would  be 
wind  bearing  with  respect  to  target,  in  the  same  manner  as  sun  direc¬ 
tion  is  defined. 

The  data  show  no  anomalies,  and  the  first  probit  run  gives 

2 

a  total  success  with  excellent  x  values, 

The  plotted  thresholds  exhibit  an  undistinguished  scatter  about 

the  normal  range,  A  fit  of  y  =  a  type  produces  a  a  6,556. 

A  t  test  shows  that  this  distribution  is  not  significantly  dlffer- 

2 

ent  from  the  normal  range  and  a  x  shows  the  points  normally 

distributed  about  the  normal  range  with  0.95  ^  p(x  )  /i  0,9^.  It  is 
necessary  to  disregard  any  effect  of  this  condition. 
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4*13  Visual  Aid  (S«e  Appendix  M) 

This  is  the  first  of  the  conditions  that  must  be  treated  dis- 
cretel7,  and  it  has  the  greatest  effect,  considering  the  spread 
between  the  high  and  low  thresholds.  Binoculars  show  a  great  advan¬ 
tage  from  a  threshold  standpoint,  but  this  must  be  a  questionable 
advantage  when  onlj  0.01  of  the  sightings  are  made  with  them.  The 
interpretation  might  wall  be  that,  if  an  observer  sees  a  target  with 
binoculars,  he  sees  it  farther  away.  Probably  a  more  factual  inter¬ 
pretation  is  that  observers  do  not  consider  binoculars  valuable  as  a 
search  aid. 

The  difference  between  no  visual  aid  and  the  use  of  simglasses 
is  very  important  however.  Many  people  feel  that,  since  sunglasses 
protect  the  eyes  from  glare  and  strong  light  and  are  more  restful, 
they  are  an  aid  to  vision.  The  advantage  of  eye  relief  is  evidently 
outweighed  by  the  attenuation  of  contrast  in  the  colored  lens  and  the 
consequent  threshold  reduction.  A  t  test  shows  a  strong  slgnifi- 
oanoe  here  and  it  is  neoessary  to  accept  the  definite  difference 
between  using  and  not  using  sunglasses.  The  results  are  tabulated? 


T  f 


Binooulars 

8.908 

1.350 

No  aid 

6.637 

1.002 

Sunglasses 

6.265 

0.949 
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4,14  Range  Determination  Method  (See  Appendix  N) 

The  reaulta  of  this  probit  anal7sls  are  somewhat  surprising. 

It  had  not  been  considered  that  the  means  by  whloh  sighting  range  was 
determined  after  the  sighting  would  have  any  effect  on  the  threshold. 
However,  the  following  explanations  may  be  worth  considering.  If  an 
observer  estimates  the  sighting  range  with  radar,  the  inference  is 
that  he  has  his  radar  activated  prior  to  sighting  and  picks  up  the 
target  on  it,  thereby  having  a  strong  clue  as  to  where  to  look.  His 
advantage  over  the  uninformed  observer  is  obvious.  The  other  expla¬ 
nation,  possibly  valid,  is  that  while  time-distance  cheeks  (calcula¬ 
ting  distance  by  speed,  and  time  to  reach  target)  are  quite  accurate, 
and  this  is  borne  out  by  the  threshold  developed,  the  unaided  obser¬ 
ver  tends  to  \indorestimate  distances,  A  tabulation  of  results 
follows 1 


T 

f 

Radar 

8.611 

1.298 

Time-distance 

6.629 

1.001 

Estimate 

6.U7 

0.923 
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4.15  Type  of  Observing  Aireraft  (See  Appendix  P) 

The  results  here  are  interesting  and  may  well  measure  the 
effeotiveness  of  the  types  of  sighting  oraft.  However  missions  and 
method  of  operation  should  be  oGnsldered,  Considering  operating 
altitude,  the  measure  of  effectiveness  becomes  weaker.  The  operating 
altitudes  are  not  normally  distributed  but  are  mors  nearly  log-normal. 
The  log  normal  means  are  computed  and  the  following  table  shows 
problt  results  with  a  consideration  of  altitude  addedi 


Type 

T 

f 

Mean 

Altitude 

T(Altitude) 

(Expected) 

T/T(Alt.) 

Patrol 

6.694 

1.014 

1278 

6.2;7 

1.072 

Utility 

5.836 

0.884 

1083 

6.107 

0.956 

Helicopter  4.836 

0.733 

606 

5.778 

0.836 

While  this  analysis  is  not  too  precise,  the  T/T(Alt.)  measure  is 
slightly  more  generous  to  the  utility  plane  emd  the  helicopter  than 
is  the  f  measure.  Again,  a  oonsideratlon  of  missions  might  alter 


these  indications. 
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4tl6  Time  of  Day  (See  Appendix  P) 

Thla  condition  is  defined  for  day,  twlll^hti  and  nighty  elnoe 
eun  altitude  takes  oare  of  lighting  variation  during  the  day.  The 
usefulneeB  of  the  thresholds  of  this  condition  are  open  to  question 
since  the  t  teat  at  the  0.05  level  shows  a  slgnlfioant  difference 
between  day  and  twilight  but  not  between  day  and  night.  The  thres¬ 
holds  for  twilight  and  night  may  be  taken  as  the  same.  However  the 
populations  sampled  may  be  different.  At  night  It  is  probable  that 
many  of  the  sightings  are  of  lights,  not  ships  as  in  daytime.  In 
view  of  the  day- twilight  significance,  the  factors  are  tentatively 
accepted.  This  subject  merits  further  study. 

The  computed  valuea  are  tabulated  as  follows: 


T 

f 

Day 

6.614 

1,002 

Night 

5.583 

0.846 

Twilight 

5.496 

0.833 
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4,17  Observer  (See  Appendix  Q) 

Here  the  pilot  and  copilot  appear  to  observe  equally  well, 
and  this  is  confirmed  by  the  t  test  with  a  p(t)  of  0.920. 

There  is  a  significant  difference  between  pilot-copilot  and  bow  look¬ 
out.  This  night  seem  unusual  unless  one  considers  the  situation  of 
the  bow-lookout  and  his  excellent  view  almost  straight  down.  It  is 
a  great  temptation  to  concentrate  on  the  area  directly  under  him  when 
ha  realizes  he  sees  better  there.  As  far  as  the  waist  lookout  is  con¬ 
cerned,  he  would  be  expected  to  see  less  well  than  the  pilot-copilot 
since  he  is  looking  to  the  side.  The  data  from  the  section  on  rela¬ 
tive  bearing  of  the  target  shows  that  at  90  degrees  he  would  expect 
to  have  a  threshold  at  about  6.2  miles.  He  is  still  seeing  less  than 
this,  though  perhaps  not  significantly  so, 

A  table  of  thresholds  aiid  factors  follows: 


T 

f 

Pilot 

6.603 

1.0006 

Copilot 

6.616 

1.002 

Waistlookout 

5.853 

0.887 

Bowlookout 

5.553 

0.841 

Ihe  low  number  of  three  sightings  by  the  tall  lookout  reinforces 
the  previously  inferred  suggestion  that  this  search  position  might 
well  be  abandoned.  Tail  observations  were  n6t  conside:-ed  here. 
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4.18  Station  (See  Appendix  R) 

The  d&ta  are  presented  he]:*eln  but  no  analysis  was  attempted. 
There  are  too  many  extraneous  Influences,  such  as  local  weather  and 
local  operating  procedures  and  policies,  to  permit  analysis  of  the 
effects  within  the  restricted  scope  of  this  study. 

There  is  certainly  a  fertile  field  here  for  cultivation  by  a 
full  scale  operations  research  study. 

One  small  scale  operations  study  was  carried  out  after  a  nximber 
of  very  odd  quantities  showed  up  in  the  cai-d  sorting.  These  were 
such  reports  as  a  bearing  of  540  degrees,  a  sighting  range  of  85 
miles,  cloud  cover  of  1.40  and  so  on.  Curiosity  suggested  pulling 
out  these  cards  to  see  what  they  might  have  in  common.  It  was 
first  found  that  they  all  come  from  the  same  station,  A  further 
investigation  showed  that  they  were  all  made  on  the  same  day, 

1  January  1956.  The  reader  is  left  to  his  own  inference. 
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4 <19  Tabular  Bbctraot  of  Threshold  Factors  (See  Appendix  S) 

The  factors  (thresholds  normalized  to  normal  range)  are 
included  in  Appendix  S  and  need  no  oonment. 
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5.0  The  Test  (See  Appendix  T) 

The  original  plan  was  to  present  the  preceding  analysis  of  the 
sighting  data  as  the  complete  study.  It  became  a  matter  of  ourloslty 
to  see  if  the  normalized  thresholds  could  be  used  to  forecast  sight¬ 
ing  thresholds.  The  results  of  this  curiosity  are  so  gratifying 
that  they  are  included  as  a  test  of  the  process  and  as  a  demonstration 
of  the  possibilities  that  further  and  more  detailed  work  may  develop. 
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5.1  Description 

The  entire  deck  of  cards  recording  air  sightinge  were  disarran¬ 
ged  by  sorting  on  the  second  digit  of  the  wind  azimuth  number,  a 
classification  that  should  result  in  no  pattern.  Then  20  cards  were 
drawn  at  approximately  equal  increments  of  distance  through  the 
length  of  the  deck.  This  procedure  should  give  a  random  assortment 
of  sightings  and  a  survey  of  the  reoanstituted  records  shows  this  to 
be  apparently  true.  (See  Appendix  T  -  3  ff.) 

The  factors  corresponding  to  the  oonditions  of  a  particular 
sighting  are  taken  from  the  factor  tables  in  Appendix  S  and  the 
product  of  all  of  these  and  the  normal  range  give  the  calculated 
threshold  for  that  sighting.  This  technique  is  followed  without 
consideration,  at  this  stage,  of  the  reported  range.  If  a  condition 
is  not  reported,  it  is  omitted,  which  is  essentially  Judging  that 
the  expeotad  value  of  this  condition  pertained  at  the  time  of  the 
sighting. 


i 

I 


52 


SIO  Reft  62.13 


5.2  The  Results  (See  Appendix  T) 

The  results  of  the  twenty  computations  are  plotted  on  the  graph 

against  the  reported  value.  A  mean  line  is  included,  as  is  the  Ideal 

line.  If  the  forecast  were  perfect,  the  computed  poli;ts  would  be 

expected  to  be  normally  distributed  about  the  ideal  line.  As  it 
2 

results,  a  x  shows  the  points  are  normally  distributed  about 

2 

the  mean  line  with  a  p(x  )  between  0.90  and  0.95. 

The  correlation  coefficient,  0.695,  shows  the  strength  of  the 
trend.  The  coefficient  is  not  expected  to  be  high  for  the  computed 
points  are  forecasted  thresholds,  and  sightings  based  on  the  condi¬ 
tions  determining  the  forecast  should  vary  normally  with  resepct  to 
the  threshold. 

These  results  suggest  strongly  that  thresholds  csin  be  effective¬ 
ly  forecasted  by  proper  consideration  of  the  conditions  affecting 
sighting  at  the  tine. 
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6.0  Conelusiona 

6.1  General  Obaervationa 

The  limited  scope  of  this  study  precludes  detailed  study  of  all 
the  inter-relations  suggested  by  the  work  on  the  data.  Some  of  the 
deficieneies  are  that  dependence  of  variables  is  treated  only  super¬ 
ficially;  trends  of  sub-classes  of  the  conditions  are  not  determined; 
complete  study  of  the  few  anomalies  occurring  in  the  computations  are 
not  possible  at  this  time.  However  it  appears  that  the  study  indi¬ 
cates  a  positive  method  of  forecasting  sightings  and  sighting  proba¬ 
bilities.  Farther,  it  would  seem  that,  until  more  conclusive  infor¬ 
mation  is  available,  the  factor  table  as  It  stands  might  well  be  a 
valuable  tool  in  forecasting. 

It  is  not  possible  to  estimate  the  total  number  of  possible 
sightings  in  any  particular  case,  since  sightings  not  made  generally 
cannot  be  reported.  In  a  standard  analysis  of  variance  procedure  the 
missed  sightings  would  lead  to  an  obviously  erroneous  result.  How¬ 
ever  one  must  assume  that  these  missing  possibilities  are  accounted 
for,  to  a  certain  extent,  in  probit  analysis  which  considers,  in 
fitting  a  curve  to  the  data,  the  character  of  the  distribution  func¬ 
tion  as  well  as  the  numerical  values  of  the  frequencies.  This 
assumption  is  further  reinforced  when  the  characteristics  of  the  log¬ 
normal  distribution  are  reviewed.  In  a  normal  distribution  the 
missed  sightings  would  result  in  a  cumulative  probability  of  less 
than  one  at  sero  range  which  infers  that  it  would  be  possible, 
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teohnlcally,  to  ooinpute  soim  eightinga  at  negative  range.  This 
situation  is  not  true  In  the  log-normal  case  for  the  oumulatlve 
probability  at  zero  range  must  be  one. 

It  Is  possible  that  further  development  of  the  problt  analysis 
method  will  allow  accounting  for  those  missed  sightings. 
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6.2  Value  of  the  Study 

This  study  la  expected  to  be  of  value  in  vision  and  vlaibillty 
reaearoh  as  a  factual  basis  for  checking  theoretical  studies. 

The  natter  herein  should  be  of  use  to  the  Coast  Guard  and 
Navy  in  refining  search  procedures  and  nethods  as  defined  in  the 
'’National  Search  and  Rescue  Manual"  (22). 

The  study  should  also  be  of  use  in  any  operational  research 
Involving  searching  and  visibility  at  sea. 
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6,3  Further  Studies  Suggested 

The  most  apparent  future  project  ia  a  similar  analysis  of  the 
sightings  from  surface  craft  and  a  comparison  with  these  results, 

There  is  a  fertile  field  for  anyone  interested  in  the  psycho¬ 
logical  study  of  the  choice  of  values  when  estimating  quantities  and 
related  fields.  In  some  processes  obsernrers  tend  to  prefer  the  num¬ 
ber  seven  to  eight,  while  in  others,  the  reverse  is  true,  It  is 
understood  that  much  research  has  been  done  on  this  subject,  and  this 
basic  data  should  yield  a  mass  of  supportive  evidence. 

Further  analysis  of  perturbations  in  the  distributions  other 

than  those  caused  by  choice  of  numbers  would  be  desirable.  The  fact 

that,  in  distributions  involving  large  numbers  of  sightings,  high 
2 

X  values  generally  occur  would  indicate  that  in  the  limit  these 
populations  might  wander  from  the  log-normal. 

Additional  investigation  would  probably  determine  asymptotes 
at  extremes  of  some  of  the  variables,  such  as  altitude  and  meteorolo¬ 
gical  visibility. 

Further  investigation  should  be  made  of  twilight  and  night 
sightings  with  perhaps  a  data  collection  program  incorporated. 

The  odd  curvature  Indicated  by  the  altitude  data  is  intriguing 
and  suggests  an  attempt  at  determination  of  whether  these  two  points 
of  inflection  are  real  or  only  coincidental  vagaries.  It  might  be 
that  there  is  a  haze  layer  or  other  condition  generally  existent 
which  causes  such  an  effect. 
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Further  study  should  be  made  of  the  few  erratlo  thresholds 
noted  in  the  disousaion  to  determine  whether  or  not  they  are  sinqple 
manifestations  of  probability  theory. 

Finally,  one  definite  reoommendationi  the  Coast  Guard  and  other 
search  agenoiea  should  give  serious  thought  toward  eliminating  tail 
lookouts,  since,  as  was  noted  above  in  Sections  4.8  and  4»17,  the 
probability  of  stern  sightings  is  close  to  zero. 
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UNITED  STATES  COAST  GUARD 


COMMANDANT 

UA.  COAAT  aUAAD 
HCADOUAimCM 
WASHINSTON  U.  D  C. 


OPmTIOHS  INSTRUCTION  NO.  58-55 

Subjt  Sighting  Data  Report  (Form  00-3627) j  instructiona  for 

1,  Purpoae,  To  prescribe  procedures  which  are  required  of  aircraft  and 
certain  floating  units  relative  to  the  preparation  and  submission  of  data 
collected  in  connection  with  the  program  for  the  collection  of  sighting 
data, 

2,  Objective.  This  program  is  designed  to  collect  reports  of  8-10,000 
sightings  of  life  rafts,  emergency  visual  signals,  small  boats  and  vessels 
under  many  visibility  and  air  and  sea  conditions, 

3,  Information.  Presently  available  "Effective  Visibility"  tables  do  not 
inoluiie  small  boats  and  vessels  with  which  the  Coast  Ouarcl  is  coimonly  coiv 
eemed,  nor  is  the  condition  of  air  and  sea  taJcen  into  consideration*  There¬ 
fore,  in  order  to  obtain  more  realistic  tables  on  this  important  subject, 
the  U,  S.  Navy,  at  the  request  of  the  Coast  Guard,  has  agreed  to  yfaluate 
(by  use  of  Univac  machines)  sighting  data  collected  by  the  Coast  Guard  and 
to  derive  empirical  formulae  from  which  curves  for  search,  sweep  width, 

and  sighting  effectiveness  may  be  drawn.  These  resxilts  will  \iltiinately  be 
incorporated  in  a  Coast  Guard  Search  and  Rescue  Manual, 

I),  Action. 

a.  Floating  units  83'  in  length  and  over  and  aircraft  shall  fill  In 
subject  form,  which  is  self  explanatory,  on  each  sighting  deemed  to  be  ad¬ 
vantageous  to  the  program.  Data  must  be  complete  for  each  sighting  reported. 
Forma  should  be  carried  on  all  flights  over  water  and  on  bridges  of  floating 
units  ready  for  use  as  may  be  practicable, 

b.  Units  shall  submit  forms  to  Commandant  (0)  in  lots  of  100  sighting 
reports, 

5,  Availability  of  Fomui.  An  initial  distribution  of  Form  CQ-3627  will  be 
made  in  the  near  future  to  all  aviation  units  and  floating  units  83'  in 
length  ar«d  ovor,  Tiw  form  will  be  included  in  ttie  Catalog  of  Forms  (CO  218) 
with  source  of  supply  "SC", 
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COPT 


OPERATIONS  INSTRUCTION  NO.  58-55 


6.  Effeotiv  date.  This  Instruction  is  effeetlTe  upon  receipt  and  uUl 
be  cilicAiba  by  separate  instruction  upon  eonpletion  of  the  project  > 


H.  C.  mXINS 
By  direction 

Enols  (1)  Sighting  Oats  Report, 

Form  CO-3627 

Diet.  (SCI,  No.  61) 

At  a,aabcd(5)}  efi(3)l  g.1.2.3,  hjklmnd) 

Bt  C(l5)j  eehi(5)!  jl(3)}  d(2)}  b(l) 

Ct  A(S)s  b4(3) 

Di  HONE 
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FLOW  OF  PROBIT  CALCUUTION 


The  notation  herein  io  defined  for  this  appendix  only. 


1.  Given  (x^,  n^,  r^)  ,  etimulue  levelsi  number  of  presentations 
and  number  of  reeiponaes  at  each  level,  determine  frequency, 

p  =  ^  .  Notet  Hereafter  the  subacripta,  1  ,  are  to  be 
Inferred. 

2.  Determine  experimental  probits,  EP,  such  thatt 


„  _  /  1  -  V2  U" 

P  -  2Tr  ® 

/  -  00 


3.  Determine  beat  fit  trial  line  and  trial  problts  I  by  least 
square  method  giving  approximate  regression  of  EP  on  x  byt 


T  s  a  +  bx 


z 

4.  Determine  weights  w  = 

Z  -  e’ 

^  2n  ® 


l/2n  e” 


■*l/2 


du;  Q  s 


5.  Determine  working  probits  y  =  P-  y  +  T 


6.  Determine  Znw,  2nwx,  2nwy,  Zhwx^^  Siwy^,  2>iwxy 
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7,  Determine 

X 

Enw 

V  = 

^  Enw 

8.  Now: 

®xx 

= 

Enwx^ 

1 

'll 

^xy 

= 

Enwxy 

Enwx  Enwy 

Enw 

®yy 

= 

Snwy^ 

.  S-SfflL. 

Enw 

9.  Then: 

‘’i 

= 

5.^ 

and  a^^  =  y  -  bx 

10.  +  bj^x.  The  subseripta  hare  and  in  9.  above 

Indicate  a  next  approxination. 


11.  The  X  applied  and,  if  not  aatlofactory,  recycling 

begins  at  3.  above  using  a^^  and  in  place  of  a  and  b  . 

12.  When  iteration  proves  satisfactory  the  following  are 
determined. 


-a 


Threshold ,  T  =  -j— 


Variance, 


Varieuice  of  threshold, 


Variance  of  a, 

ft  O 

XX 

Variance  of  b,  S?  = 

^xx 


Enw 


*  M 


XX 


Variance  of  standard  deviation,  S 


=2  - 


XX 
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NOTATION  j 

i 

i 

A 

Altitude,  1000*8  feet  ■ 

1 

j 

B 

Target  bearing,  degrees 

C 

Cloud  cover,  decimal  fraction 

! 

i 

D(P) 

Noi*mal  deviate  for  probability  p 

f 

Factor 

F(o) 

Factor  function  of  .  ,  any  variable 

L 

Length  of  vessel,  feet 

OB 

Obsejn^er 

OU 

Observing  unit 

P 

emulative  probability  of  sighting 

P(«) 

Probability  of  .  ,  any  function 

r 

Correlation  coefficient 

RD 

Range  determination  method 

S 

Height  of  major  swells,  feet 

B 

Standard  deviation 

i 

Mean  standard  deviation 

®8 

Standard  deviation  of  8 

Standard  deviation  of  S 

Standard  deviation  of  T 

SA 

Sun  altitude,  degrees 

SB 

Sun  bearing,  relative  to  observer-target  line,  degrees 

ST 

Air  station 

T 

Threshold,  miles 

1 

1 

Tn 

Normal  range,  miles  APPENDIX  X 

SIO  Ifleft  62-13 


lao 

T(*)  Threshold  foaotion  of  .  ,  any  variable,  tallea 

TD  Time  of  day,  general 

t  Student's  t  funotiwi 

V  Vlalblllty,  miles 

VA  Visual  aid 

WA  Wind  azimuth,  degrees 

WS  Wake  size,  proportional  to  length  of  vessel 

WV  Wind  velocity,  knots 

2 

X  The  accumulated  value  of  the  measure  of  goodness 

of  fit  and  other  statistical  tests. 
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NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  puirose 
other  than  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rl^ts 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 
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BALLISTIC  RESEARCH  LABORATORIES 


Id 

ABERDEEN  PROVING  GROUND,  MARYLAND 


ASTIA  AVAILABILITY  MOTT.CE 


Qualified  requestora  may  obtain  copies  of  this  report  from  ASTIA. 


The  findlnga  in  thj.s  report  are  not  to  be  construed 
as  an  official  Department  of  the  Army  position. 
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THE  SIMULATION  OR  INTERIOR  BALLISTIC  PERFORMANCE  OF 
GUNS  BY  DIGITAL  COMPUTE  PROGRAM 

Paul  G.  Baer 
Jerome  M.  Frankie 

Interior  Ballistics  Laboratory 


RDT  8s  E  Project  No.  1M010501A004 

ABERDEEN  PROVING  GROUND,  MARYLAND 


BALLISTIC  RESEARCH  LABORATORIES 

REPORT  NO.  1183 


PGBaer/JMPrankle/mec 
Aberdeen  Proving  Ground,  Mi. 
December  1962 


THE  SIMULATION  OF  INTERIOR  BALLISTIC  PERFORMANCE  OF 
GUNS  BY  DIGITAL  COMPUTEK  PROGRAM 

ABSTRACT 

When  non- conventional  guns  are  to  be  considered  or  when  detailed  design 
Infoimtion  is  required,  interior  ballistic  calculations  become  more  difficult 
and  time-consuming.  To  deal  with  these  problems,  the  equations  which  describe 
the  interior  ballistic  performance  of  giuis  and  gun-like  weapons  have  been  pro*' 
grammed  for  the  high-speed  digital  computers  available  at  the  Ballistic  Research 
Laboratories.  The  major  innovation  contained  in  the  equations  derived  in  this 
report  is  the  provision  for  use  of  propellant  charges  made  ut  of  several  pro¬ 
pellants  of  different  chemical  compositions  and  different  grsnulatlone .  Re- 
svilts  obtained  by  the  method  described  in  this  report  compare  favorably  with 
those  of  other  interior  ballistic  systems.  In  addition,  considerably  more  de¬ 
tail  is  obtained  in  far  less  time.  A  comparison  with  experimental  data  from 
well- instrumented  gw-flrings  is  also  presented  to  demonstrate  the  validity  of 
this  method  of  computation. 
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a  acceleration  of  projectile.  In. /sec’ 

a  constant  defined  by  Equation  (28a),  dimensionless 
0 

A  area  of  base  of  projectile  Including  appropriate  portion  of  rotating 
band,  in.^ 

■z 

bj^  covwj-ume  of  1  th  propellant,  in./Xh 
c  diameter  of  bore,  in. 

specific  heat  at  constant  volume  of  1  th  propellant  (c^  la  a  function 
^  of  T),  in.-lb/lb-°K  ^ 

—  mean  value  of  specific  heat  at  constant  volume  of  1  th  propellant 
^1  (over  temperature  range  T  to  T^  ),  in.-ib/lb-*^K 

—  mean  value  of  specific  heat  at  constant  pressure  of  1  th  propellant 

C  Q 

^1  (over  temperature  range  T  to  T^  ),  In.-lb/lb-  K 
Initial  weight  of  1  th  propellant,  lb 
Cj.  Initial  weight  of  Igniter,  lb 

d^  diameter  of  perfoi-ation  in  1  th  propellant  grains,  In. 
dt  Incremental  time,  sec 

o 

dT  incremental  temperature,  K 

dx  incremental  distance  traveled  by  projectile,  in. 

(1.Z  "I 

1  mass  fraction  burning  rate  for  1  th  propellant,  sec 
dt 

Dj,  outside  diameter  of  1  th  propellant  grains,  In. 

cnorsi'  Icct  due  to  heat  looc.  In. -lb 

E  kinetic  energy  of  propellant  gas  and  unbumod  propellant,  in. -lb 
P 

Epy  energy  lost  due  to  bore  friction  and  engraving  of  rotating  band,  In, -lb 
functional  relationship  between  and 
F^  resultant  axial  force  on  projectile,  lb 


7 


Fj  frictional  force  on  projectile,  It 

F^  "force"  of  i  th  propellant,  in.-lt/lt 

Fj  "force"  of  igniter  propellant,  in.-lb/lt 

F  propulsive  force  on  base  of  projectile,  lb 
P 

F  gas  retardation  force,  lb 
r 

2 

g  constant  for  conversion  of  veight  units  to  mass  units,  in. /sec 

G  functional  relationship  between  p^  and  x 

K  burning  rate  velocity  coefficient,  in, 

^  sec  in. /sec 

K  burning  rate  displacement  coefficient,  in. 

^  sec- in. 

length  of  1  th  propellant  grains,  in. 

specific  mass  of  i  th  propellant,  Ib-mcls/mol 

M  mass  of  projectile,  alugs/lS 

n  number  of  propellants,  dimensionless 

n*  ratio  defined  by  Equation  (28b),  dimensionless 

munber  of  perforations  in  1  th  propellant  grains,  dimensionless 

p  space-mean  pressure  resulting  from  burning  i  propellants,  psi 

pressure  on  base  of  projectile,  psi 

p  pressure  of  gas  or  air  ahead  of  projectile,  psi 
S 

space-mean  pressure  resulting  from  burning  of  1  th  propellant,  psi 
p^  igniter  pressure,  pal 
breech  pressure,  psi 
p^  resistance  pressure,  psi 

Q  .  energy  released  by  burning  propellant,  in. -lb 
linear  burning  rate  of  i  th  propellant,  in. /sec 

adjusted  linear  burning  rate  of  1  th  propellant,  in. /sec 
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functional  relationship  between  and  p 

2 

surface  area  of  partially  burned  i  th  propells»t  graift,  in» 

2 

surface  area  of  an  unburned  i  th  propellant  grain,  in. 
time ,  sec 

0 

mean  temperature  of  propellant  gases,  K 
adiabatic  flame  temperature  of  i  th  propellant,  °K 

adiabatic  flame  temperature  of  Igniter  propellant,  °K 

temperature  of  unburned  solid  propellant,  °K 

two  times  the  distance  each  surface  of  i  th  propellant  grains  has  receded 
at  a,  given  time,  in. 

Internal  energy  of  propellant  gases,  in. -lb 
velocity  of  projectile,  in. /sec 
velocity  of  projectile  at  muzzle  of  gun,  in. /sec 
specific  volume  of  propellant  gas,  In.^/lb 

3 

volume  behind  projectile  available  for  propellant  gas,  in. 

3 

volume  of  an  unbxirned  i  th  propellant  grain,  in. 

3 

volume  of  empty  gun  chamber,  in. 
external  work  done  on  projectile,  in. -lb 
weight  of  projectile,  lb 
travel  of  projectile,  in. 

travel  of  projectile  when  base  reaches  muzzle,  in. 


fraction  of  mass  oi 


^  4  4'V« 


propellant  bur.ned,  dimensionless 


fraction  of  mass  of  igniter  burned,  dimensionless 
bvurning  rate  exponent  for  i  th  propellant,  dimensionless 


burning  rate  coefficient  for  1  th  propellant, 


in.  1 
■  psl“ 


effective  ratio  of  specific  heats  as  defined  by  Equation  (27a),  dloehsioaleas 
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ratio  of  specific  heats  for  i  th  propellant,  dlmenBionless 
ratio  of  specific  heats  for  igniter  propellant,  dimensionless 
Pidduck-Kent  constant,  dimensionless 
density  of  i  th  solid  propellant,  Ib/in."* 
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ITOODUCTIOK 


The  interior  'balllstlcian  must  frecLuently  predict  the  interior  ballistic 

perfomnance  of  guns.  In  some  instances,  it  is  sufficient  to  calculate  muzzle 

velocity  and  maximum  chamber  pressure  for  a  conventional  gun  from  a  knowledge 

of  the  propellant  charge,  the  projectile  wei^t,  and  the  gun  characteristics. 

f  1  w 

This  calculation  is  usually  referred  to  as  the  classical  central  problem  '  ' 
of  interior  ballistics.  When  non-conventional  guns  are  considered  or  when 
detailed  design  Information  is  required,  it  is  necessary  to  know  more  than 
these  two  salient  values.  For  the  more  demanding  problems,  complete  interior 
ballistic  trajectories  may  have  to  be  calculated.  These  trajectories  consist 
of  displacement,  velocity,  and  acceleration  of  the  projectile  and  chamber  pres¬ 
sure,  all  as  functions  of  time. 

The  literature  of  interior  ballistics  contains  descriptions  of  many  methods 
for  solving  the  problem  of  predicting  the  performance  of  guns.  Methods, 

varying  from  the  piurely  en5)irical  to  the  "exact"  theoretical,  have  been  devised 
in  tables,  graphs,  nomograms,  slide  rules,  and  sinqjlified  equations  solved  in 
closed-form.  Some  of  these  methods  require  data  from  the  firing  of  the  gun 
being  considered  or  from  a  very  similar  gun.  All  of  these  methods  require  some 
simplification  of  the  basic  equations  of  interior  ballistics. 

To  eliminate  the  restrictions  imposed  by  assuniptionB  made  only  to  facil¬ 
itate  the  mathematical  solution  of  the  problem,  the  interior  ballistic  equations 
have  been  programmed  for  high-speed  electronic  computers.  Both  analog  and  dig¬ 
ital  oongjuters  have  been  used  to  calculate  detailed  interior  ballistic  trajec¬ 
tories.  There  are  advantages  and  disadvantages  associated  with  each  type  of 

f  3) 

computer.  Several  years  ago,  '  the  interior  ballistic  equations  were  pro¬ 
grammed  for  the  digital  computers^**  available  here  at  the  Ballistic  Research 
Laboratories.  Since  that  time,  considerable  use  has  been  made  of  this  program 
for  studying  gun  and  gun-llke  systems  and  for  routine  calculations. 


*  Superscripts  indicate  references  listed  at  the  end  of  this  report. 

**  Althougjh  the  interior  ballistic  equations  were  originally  programmed  only 

(4) 

for  the  ORDVAC,  '  '  they  have  been  recently  reprogrammed  In  more  general 
form  for  the  ORDVAC  and  the  newer  BRLESC. 


1!he  computer  program  descrl'bed  in  this  report  has  been  designed  tg  solve 
a  set  of  non-linear,  ordinary  differential  and  algebraic  equations  which  simu¬ 
late  the  Interior  ballistic  performance  of  a  gun.  In  this  method,  the  usual 
set  of  equations  which  pertains  to  the  burning  of  a. single  propellant  has  been 
modified  to  account  for  the  burning  of  composite  charges,  l.e.,  charges  made 

up  of  several  propellants  of  different  chemical  compositions  and  different 
* 

granulations.  The  computer  program  m^  be  suitably  modified  to  study  non- 
conventlonal  guns  and  gun-like  systems.  A  number  of  these  optional  programs 
have  been  devised  and  used  extensively. 

INTERIOR  BALLISTIC  THEBRY 

Interior  Ballistic  System 

The  basic  components  of  the  interior  ballistic  system  for  a  conventional 
gim  are  shown  in  Figure  1.  A  set  of  equations  can  be  formulated  which  math¬ 
ematically  describes  the  distribution  of  enerQr  originating  from  the  burning 


Figure  1.  Basic  Components  of  the  Interior  Ballistic 
System  for  a  Conventional  Gun 


*  The  present  program  can  be  operated  with  as  many  as  five  different  types  of 
propellant  charges  for  each  problem. 

**  See  Section  entitled  Options  to  Routine. 
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propellant  and  the  suhseqLuent  motion  of  the  components  of  the  system,  in  the 
development  which  followsi  two  major  assumptions  are  made  to  accowt  for  the 
behavior  of  composite  charges: 

1.  The  total  chemical  energy  available  is  the  simple  sum  of  the  chemical 
energies  of  the  Individual  propellants. 

2.  The  total  gas  pressure  is  the  simple  sima  of  the  "partial  pressures" 
resulting  from  the  burning  of  the  individual  propellants. 

Energy  Equation 

Application  of  the  law  of  conservation  of  energy  leads  to  the  energy  equa¬ 
tion  of  interior  ballistics.  This  may  be  written  as: 

Energy  Released  _  Internal  Energy  ^  External  ^  Secondary 

by  Burning  Propellant  of  Propellant  Oases  Work  Done  Energy  Losses 

on  Projectile  (l) 


or: 


Q  =  U  +  W  +  Losses  (la) 

In  Equation  (la)  the  energy  released  by  the  burning  propellant  (q)  Is  assumed 
to  be  equal  to  the  simple  sum  of  the  energies  released  by  the  individual  pro¬ 
pellants  as  previously  stated.  Therefore: 

1=1  ^  (2) 


Because  of  gas  expansion  and  external  work  performed  In  a  gun>  the  gas  temper¬ 
ature  is  less  than  the  adiabatic  flame  temperature  (T^  ).  The  Internal  energy 
of  the  gas  (U)  is  then:  ^ 


n 


ial 


The  external  work  done  on  the  projectile  is  given  by: 


X 

W  =  A  / 
*^0 


Pb  <3^ 


(3) 


(4) 
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Substituting  Equations  (2),  (3),  and.  (4)  into  Equation  (la)  gives: 

Z  [°l"l  I  =  Z  [°l*i  ^  i 


which  may  be  rewritten  as; 


f  ^  c  dll  -  A  f  p,  dx  +  Losses 
±  ±  }  Jq 

As  the  c  do  not  vary  greatly  over  the  temperature  ranges  from  T  to  T  , 
^i  °i 

they  can  be  replaced  with  mean  values  (o  ).  Integration  of  Equation  (5) 

.  '^i 

gives : 


n  ^ 

ZC.z  0  (T  -  T)  =  A  /  p, 
i  1  Vf  ^  o^  Jq  ^b 


dx  +  Losses 


and  solving  for  T: 


“  3(. 

Z°A  V\  -  ‘'i  h 

i=l _ ^  ^ _ 

n 

Z  \ 

J  1  1 


dx  -  Losses 


Next,  the  "force"  of  each  propellant  is  defined  by: 


F,  =  ra,RT 
i  i  0. 


and  the  well-known  relations: 


cpi  -  ovj^  = 


are  Introduced. 


Combination  of  Equations  (9)  and  (lO)  gives: 


(7^  -  1) 


m.R 

1 


Substitution  of  Equation  (ll)  into  Equation  (8)  gives: 

F. 


°1  (7j  -  1)  \ 


■1  (12) 

Finally,  substitution  of  Equation  (l2)  into  Equation  (7)  gives  Rfesal's  equation 
in  the  form: 


F,C,2 


dx  -  Losses 


T  » 


1=1 


F.C  -2, 


V  ill 

L  TypTrX 

1=1  ^  ° 


(13) 


For  moat  problems,  it  is  convenient  to  assume  the  igniter  completely  burned 
(zj.  =  1)  at  zero- time.  Equation  (l5)  may  be  restated  as: 

F  C  z  1  F  C  ^ 

-  + - ~  ^  I  ^b 

-  ^  0 

y 

L  1  =  1  <'1 


r  n 


.Vl 


(14) 


The  terms  ^  J  ^  Losses  of  Equation  (l4)  can  now  be  considered 

in  more  detail.  The  work  done  on  the  projectile  results  in  an  equivalent  gain 
in  kinetic  energy  of  the  projectile  except  for  losses.  Including  these  losses 
under  the  general  category  of  energy  losses: 


W. 


^  -I 


(15) 
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l2) 

Accordiog  to  Hunt,  '  '  the  energy  lossee  to  be  considered  are: 

(1)  kinetic  eneror  of  propellant  gas  and  unburned  propellant, 

(2)  kinetic  energy  of  recoiling  parts  of  gun  and  carriage, 

(3)  heat  energy  lost  to  the  gun, 

(4)  strain  energy  of  the  gun, 

(5)  enerQT  lost  in  engraving  the  rotating  band  and  in  overcoming  friction 
down  the  bore, 

and 

(6)  rotational  energy  of  the  projectile. 

For  discussion  of  each  type  of  secondary  energy  loss,  see  Reference  (2). 

Types  (2),  (4),  and  (6)  are  estimated  to  be  less  than  one  percent  for  each 
category  and  have  been  neglected  here. 

The  kinetic  energy  of  propellant  gas  and  unburned  propellant  can  be  repre- 
sented  by  '  ^ 

(i ..)  - 

P  ~  2g6  (16) 

The  energy  losses  resulting  from  heat  lost  to  the  gun  can  be  estimated  by  a 

(2) 

semi -empirical  relationship  described  by  Hunt:' 


n 


At  the  present  time,  the  Introduction  of  a  more  sophisticated  treatment  of  heat 
loss,  with  its  attendant  complexity,  does  not  seem  to  be  warranted.  Such  a 
substitution  can  be  made  if  and  when  It  appears  desirable. 


The  final  energr  losses  to  be  considered  here  consist  of  those  resulting 
from  engraving  of  the  rotating  band,  friction  between  the  moving  projectile 
and  the  gun  tube,  and  acceleration  of  air  ahead  of  the  projectile.  Individual 
estimates  of  these  are  difficult  to  make,  so  they  have  been  grouped  as  resis¬ 
tive  pressure  in  the  form: 


(18) 


The  versus  x  function  is  discussed  In  greater  detail  in  the  section  concern¬ 
ing  forces  acting  on  the  projectile. 


Substitution  of  Equations  (15)»  (l6)»  (17)>  and  (iQ)  into  Equation  (iJj-) 
results  in  the  form  of  the  energy  equation  used  in  this  computer  program: 


Equation  of  State 


The  pressure  acting  on  the  base  of  the  projectile  cem  be  calculated  from 
a  series  of  equations,  once  the  temperature  of  the  gas  is  determined  from  the 
energy  equation.  Generally,  the  equation  of  state  for  an  ideal  gas  takes  the 
form: 


^l\  “  “l 

where  =»  the  volume  per  unit  mass  of  1  th  propellant  gas. 

Mow,  define  the  volume  behind  the  projectile  which  is  available  for  pro¬ 
pellant  gas,  as: 
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.  rj!  S-j:KkiW(#1I 


Volume  Available 
for  Propellant 
Gas 


Initial  Empty 
Chamber  Voliame 


Volume  Occupied 
by  Unburned 
Solid  Propellant 


Volume  Resulting 
from  Projectile 
Motion 

Volume  Occupied 
by  Gas  Molecules 
(covolume)  /qt 


VaV+Ax-Vc.  /-  vV 
CO  L,  (1-zJ  -il/  ill 


l=>l  P 


By  the  definitions  of  Equations  (20)  and  (21), 


Substituting  Equations  (8)  and  (23)  into  Equation  (20)  and  rearranging  gives: 


F.C.z.  T 

®1  V  T 
c  o^ 


If  the  b^  are  assumed  to  be  constants  over  the  tenrperature  range  from  T  to 
,  and  if  the  total  gas  pressure  is  taken  as  the  simple  sum  of  the  "partial 
prlssures"  resulting  from  the  burning  of  the  Individual  propellants  as  pre¬ 
viously  stated,  then: 


j  .  *1  » _  * 


c  1=1 


As  before,  if  it  is  assumed  that  the  igniter  is  completely  burned  (z^.  =  l) 
at  zero-time,  Equation  (25)  may  be  restated  as; 


Wi.)  * 

°  L'&i  Oi  '  Op  J 


1 


The  space-mean  pressure,  p,  given  ty  Equation  (26)  is  used  in  the  cal¬ 
culation  of  the  fraction  of  propellant  hurned  at  any  time.  This  relationship 
is  discussed  in  the  section  concerning  hurning  rates.  There  Is,  however,  a 
pressure  gradient  from  the  breech  of  the  gun  to  the  base  of  the  projectile 
which  must  be  considered  in  developing  the  equations  of  motion  for  the  pro¬ 
jectile.  This  pressure-gradient  problem  was  first  considered  by  Lagrange  and 
is  commonly  referred  to  as  the  Lagrange  Ballistic  Problem.  Later  studies  in 
this  area  were  made  by  Love  and  Pidduck,  Kent,  and  others.  For  this 
computer  program,  the  improved  Pidduck-Kent  solution  developed  by  Vlnti  and 
Kravltz  has  been  used: 


P 


.1=1 

1  +  - 


V  (27)* 

In  addition  the  breech  pressure,  p^,  is  calculated  by  the  method  contained 
in  Reference  (6).  This  is  the  pressure  usually  measured  in  experimental  inte¬ 
rior  ballistic  studies: 


where :  l/a 

o 


2  n’+5  2  (n’+l) 


i=l 


(28) 

(28a) 


In  Reference  (6),  t.he  detenral natl nn  of  B  depends  on  the  ratio  of  epeftiflc 

heats,  7.  For  composite  charges,  an  effective  value  n 

is  used  for  this  purpose.  V  _ 

^  b.7. 
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(276) 


(26b) 


1 

and  n*  =  - 

7'-l 

Mass -Fraction  Burning  Rate  Equation 

Both  the  energy  equation  (Equation  (19))  and  the  equation  of  state  (Equa¬ 
tion  (26))  are  algebraic  equations  whose  solutions  depend  upon  the  solutions 
of  several  non-linear,  ordinary  differential  equations.  The  mass-fraction  burn¬ 
ing  rate  equation  expresses  the  rate  of  consumption  of  solid  propellant  and 
hence  the  rate  of  evoluliun  of  propellant  gas.  This  may  be  written  as: 


dz 


1  = 


where : 


dt 


r,  =  R 


1 

V 


Sf  r^ 


1 


(29) 

(30) 


and;  =  f^  (s^)  (5I) 

For  most  gun  propellants,  Equation  (50)  may  be  quite  satisfactorily  stated  as: 
ri  =  Pj(p)“i 


For  certain  propellants.  Including  those  plateau  and  mesa  types  used  In 
solid-fuel  rockets.  Equation  (52)  will  not  suffice  for  gun  calculations.  In 
these  cases,  it  is  preferable  to  make  use  of  a  tabular  listing  of  r^'s  and 
corresponding  p's  (Equation  (50))  and  to  interpolate  for  the  desired  r^.  The 
r^'s  calculated  by  either  Equation  (50)  or  Equation  (52)  are  closed  chamber 
burning  rates.  As  discussed  In  later  sections  of  this  report,  these  burning 
rates  may  be  increased  by  addition  of  factors  proportional  to  the  velocity  and 
displacement  of  the  projectile  in  the  following  manner; 

=  ^i  +  Kv  V  +  X  (32a) 
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Similarly,  the  form  fvmctlon  described  by  Equation  (jl)  may  be  stated 
In  one  of  several  ways.  In  many  Interior  ballistic  systems,  the  form  function 
is  chosen  for  convenience  of  analytical  solution.  Where  routine  numerical 
con5)utations  are  handled  by  use  of  a  high-speed  digital  con^uter,  the  geomet¬ 
rical  form  of  the  propellant  grain  may  be  used  to  obtain  the  functional  rela¬ 
tionship,  fj^,  between  S^  and  z^.  For  the  usual  grain  shapes  encountered,  these 
equations  are  given  in  Appendix  A.  This  Appendix  also  contains  the  method  for 
handling  such  equations  in  the  conqjuter  routine.  To  extend  these  equations  to 
include  propellant  slivering  see  Reference  (9) • 

Equations  of  Pro.ieetlle  Motion 

The  translational  motion  of  the  projectile  down  the  gun  tube  may  be  cal¬ 
culated  from  the  forces  acting  on  the  projectile.  Figure  2  shows  the  axial 
forces  considered  in  determining  the  resultant  force. 


Figure  2.  Axial  Forces  Acting  on  Projectile 

The  propulsive  force,  P^,  is  that  resulting  from  the  pressure  of  the 
propellant  gas  on  the  base  of  the  projectile  according  to: 

=  PfeA  (55) 

where  Pj^  is  obtained  from  Equation  (27). 

The  frictional  force,  F^,  is  the  retarding  force  developed  by  resistance 
between  the  bearing  surfaces  of  the  projectile  and  the  inside  of  the  gun  tube. 
This  is  usually  the  resistance  between  the  rotating  band  and  the  rifling  of 
the  tube  and  includes  the  force  required  to  engrave  the  rotating  band.  It  may 
be  e:q)reseed  as: 

P-  “  p  A  (3^) 

r  T 
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The  determination  of  is  difficult  in  most  cases.  Many  interior 
ballistic  solutions  use  an  increased  projectile  mass  (approximately  %)  to 
account  for  its  effect.  There  are  several  disadvantages  inherent  in  such  a 
treatment.  Althou^  the  muzzle  velocity  may  be  calculated  reasonably  well, 
the  detailed'  trajectory  will  be  altered  considerably.  It  is  not  possible  to 
simulate  the  case  where  a  projectile  lodges  in  the  bore  (see  Reference  (lO) 
for  experimental  trajectories  for  this  condition).  For  this  computer  program, 
experimental  data  of  the  type  given  in  Reference  (ll)  may  be  used  by  insert¬ 
ing  a  tabulation  of  the  function: 

=  G(x)  (Jlta) 

The  gas  retardation  force,  F^,  is  that  which  results  from  the  pressure 
of  air  or  gas  ahead  of  the  projectile,  stated  as: 

F  =  p  A  (55) 

r 

where  p  is  small  enou^  to  be  neglected  except  for  very  high  velocity  systems, 

S 

li^t  gas  guns,  and  other  special  applications.  In  the  discussion  of  the  Energy 

Equation  in  the  Interior  Ballistic  Theory  Section,  p  was  considered  a  part  of  p  . 

B  ^ 

The  resultant  force  in  the  axial  direction  is  then: 

•  ■'p  -  '■f  -  '■r 


or: 


F  = 
a 


A(p^ 


(37) 


is: 


or; 


The  acceleration  of  the  projectile,  by  Newton's  second  law  of  motion, 


-  Pj 


'r> 


a 


M 

Ag  (P^  -  Pg  -  P^) 


W 


P 


(38) 

(39) 
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Since  a 


a  dt 


the  velocity  of  the  projectile  is  given  bys 


(40) 


and  the  displacement  of  the  projectile  is  given  by: 


(41) 


Summary  of  Interior  Ballistic  Equations 

The  equations  which  arc  uced  in  the  computer  program  are  now  summarized 
for  ease  of  reference. 

Energy  Equation 
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Equatloa  of  State 
r"  n 


P  =  V. 


(I^)  * 


L  1=1  “i 


FjCi 


"I  J 


(26) 


a  C. 


where;  +  Ax  -  J  —  (l-z^)  -  J 


1=1  ^1 


i=l 


(22) 
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Mass-Fraction  Burning-Rate  Equations 


®1  "i 


(29) 


^1 


(32) 


or: 
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r .  +  K  V  +  K  X 
1  V 


(32a) 
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Eauatlons  of  Pro.jecblle  Motion 


a  =  Ag  (p.  -  p  -  p  ) 
b  g  r 

W 


(39) 


m 


(41) 


COMPUTATION  ROUTINE 

The  set  of  non-linear,  ordinary  differential  and  algebraic  equations, 
summarized  at  the  end  of  the  previous  section,  simulates  the  interior  ballistic 
performance  of  a  gun  or  gun-like  system.  A  numerical  computation  routine  has 
been  devised  for  the  simultaneous  solution  of  these  equations.  The  general¬ 
ized  flow-diagram  for  the  routine  is  presented  in  Appendix  B.  Using  the  FORAST 
language,  '  the  solution  has  been  programmed  for  the  ORDVAO  and  BRLESC  com¬ 
puters. 

Preliminary  Routine 

To  reduce  computation  time  and  conserve  memory  space,  a  prellmlnaiy  routine 
has  been  introduced.  Here  all  data  required  for  the  computation  are  read  into 
the  conq>uter,  constant  groupings  (e.g.. 


P^Ci 


^i°l 


(ri-1) 


etc.,  are  calculated  and  stored 


for  subsequent  use,  and  data  to  permanently  Identify  the  computer  run  ere  print¬ 
ed  out.  A  coaplete  listing  of  required  input  data  may  be  found  in  Appendix  C. 
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\  Main  Routine 

The  min  conrputatlonel  routine  Is  presented  In  the  generalized  tXavr- 
diagram  of  Appendix  B.  To  follow  the  procedure,  consider  the  three  sequential 
phases  of  the  problem: 

Phase  I  -  Prom  time  of  Ignition  until  the  projectile  starts  to  move. 

Phase  II  -  Prom  time  of  initial  projectile  motion  until  all  propellants 
are  consumed. 

Phase  III  -  From  time  of  propellant  burnout  until  projectile  leaves  the 
gun  muzzle. 

At  the  time  of  ignition  (Phase  I  begins),  it  Is  assumed  that  the  igniter 
is  completely  burned  (zj  =  l)  and  none  of  the  other  propellants  have  started 
to  bum  (all  =  O).  The  space-mean  pressure,  consisting  only  of  the  igniter 
pressure,  is  calculated  from: 


Equation  (42)  la  derived  from  Equations  (19)  and  (26)  by  means  of  the  slsipli- 
fylng  ignition  assumptions  stated  above. 

The  linear  burning  rate  for  each  propellant  can  now  be  determined  from 
either  Equation  (50)  or  Equation  (32)  in  combination  with  Equation  (3Sa).  If 
the  Interpolation  indicated  by  use  of  Equation  (30)  is  selected,  the  general¬ 
ized  interpolation  sub-routine*  is  employed.  Tlie  mass -fractions  burned,  z^*b, 
during  a  small  time  Interval,  dt,  are  determined  by  Integration  of  Equation 
(29).  The  surface  areas  of  the  unburned  propellant  (see  Appendix  A)  are  used 

in  this  initial  calculation*  The  Runga-Kutta  method  of  numerical  Integration, 

(12) 

as  modified  by  Glli,  '  '  is  commonly  used  for  the  solution  of  sets  of  or¬ 

dinary  differential  equations  and  has  been  employed  here. 

Calculation  of  the  tenq>erature,  T,  from  Equation  (19)  and  the  volume 
available  for  propellant  gas,  V  ,  from  Equation  (22),  will  allow  the  calcu¬ 
lation  of  the  new  space-mean  pressure,  p,  at  time,  dt,  from  Equation  (26).  The 
surface  areas  of  the  now  partially  burned  propellants  are  computed  from  equa¬ 
tions  presented  in  Appoudix  A.  All  results  of  interest  are  prlntod-out  at  this 
time  **  and  these  results  used  as  initial  conditions  for  oalculationa  during 


I 


*  See  Reference  (I6)  for  interpolation  by  divided  differences. 
**  See  Appendix  C  for  listing  of  output  data. 
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the  ensuing  time -interval.  Those  terms  in  Eauations  (17),  (19),  and  (22)  which 
Involve  velocity  or  displacement  are  zero  during  this  phase  of  the  coraputatlonf 
This  calculation-loop  is  continued  until  the  space-mean  pressure  exceeds  a 
pre-selected  "shot-start"  pressure  and  the  projectile  starts  to  move.  Phase  I, 
which  has  heen  arbitrarily  defined,  ends  at  this  time. 


Phase  II  requires  the  addition  of  the  equations  of  motion  to  the  sequence 
followed  during  Phase  I.  Equations  (27),  (59)>  (4o),  and  (4l)  are  used  to  cal¬ 
culate  the  values  of  the  acceleration,  velocity,  and  displacement  of  the  pro¬ 
jectile  at  the  end  of  each  time  interval.  Integration  specified  in  Equations 
(4o)  and.  (4l)  is  again  performed  by  the  Runga-Kutta-GlU  method.  Values  of 
velocity  and  displacement  are  now  available  for  usg  in  terms  of  Equations  (17)1 
(19),  and  (22).  To  compute  values  for  =>  A  ^  which  is  one  of 

the  terms  in  Equation  (19) >  the  generalized  Interpolation  sub-routine  must  be 
used  to  obtain  p^  from  the  tabular  Information  described  by  Equation  (54a). 
nils  integration  is  performed  by  use  of  the  Trapezoidal  Rule.* 


As  time  is  Increased  by  the  addition  of  small  time-intervals,  calculations 
during  Phase  II  are  continued  around  this  expanded  loop  with  print-out  of  appro¬ 
priate  results  at  the  end  of  each  time  interval.  One  at  a  time,  the  propellants 
are  completely  consumed  and  this  phase  is  ended.  A  series  of-  switches  has  been 
incorporated  in  the  program  to  circumvent  the  necessity  of  introducing  propel¬ 
lants  in  any  special  order.  In  fact,  it  may  not  always  be  possible  to  predict 
the  exact  order  in  which  several  different  propellants  will  be  burned  out.  The 
combination  of  the  propellant  switches  and  the  Start-of -motion  switch  makes  It 
possible  to  handle  problems  where  one  or  xoore  propellants  bum  out  before  the 
projectile  starts  to  move. 

With  all  propellants  consumed.  Phase  III  begins.  The  mass -fractions  burned 
have  all  become  unity  and  the  equations  concerned  with  burning  (Equations  (S9)> 
(51),  Btid  (50)  or  (58))  are  eliminated  from  the  loop.  As  in  the  other  phases, 


*  Although  the  Thrapezoldal  Rule  is  a  relatively  crude  method  for  numerical 
lntegr.atlon,  the  accuracy  of  the  p^  versus  x  data  available  does  not  warrant 

a  more  accurate  and  hence  more  complex  method. 
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results  are  prlnted-out  at  the  end  of  each  time-interval.  A  continual  check 
is  made  of  the  displacement  of  the  projectile  to  detcnnine  whether  or  not  it 
has  reached  the  muzzle  of  the  gun.  When  the  projectile  passes  the  muzzle, 

Phase  III  has  ended  and  the  program  is  stopped. 

It  is  possible  for  the  projectile  to  reach  the  muzzle  (and  the  program 
stopped)  before  Phase  II  is  completed.  This  would  simulate  a  gun-firing  in 
which  unburncd  propellant  is  ejected  from  the  muzzle.  It  is  also  possible 
for  the  pi-ogL'am  to  simulate  a  firing  in  which  the  projectile  becomes  lodged 
in  the  tube.  In  this  case.  Phase  HI  Is  not  completed  and  the  program  is 
stopped  wlicn  the  projectile  displacement  docs  not  increase. 

At  each  time-interval  after  the  beginning  of  Phase  II,  the  breech  pressure 
is  detemined  from  Equation  (P8)  and  printed  out.  This  result  is  not  used  in 
l.he  comiJutatlonaL  routine  but  is  used  to  compare  theoretical  and  experimental 
I’l -suits.  A  coutlnuai  check  is  made  of  the  calculated  pressures  and  the  maximum 
!u'i"ch  pressun.  Is  stored  with  Its  associated  time  and  projectile  displacement. 
This  InfoniiuLlon  is  prlnted-out  at  the  end  of  the  program.  Calculations  dur¬ 
ing  the  last  time-interval  result  In  a  projectile  displacement  somewhat  greater 
llum  the  deoii'cd  distance  to  the  muzzle.  A  linear  interpolation  between  re¬ 
sults  at  the  last  two  time-intervals  is  used  to  obtain  results  exactly  at  the 
muzzle.  These  results  are  also  prJntcd-out  at  the  end  of  the  program. 

Options  to  Routine 

A  considerable  number  of  options  have  been  designed  and  coded  for  special 
problems.  These  include  changes  which  enable  the  program  to  be  used  for  guns, 
or  gun-like  weapons,  which  are  not  of  conventional  design  (Figure  l)  and  changes 
which  vary  the  treatment  of  some  of  the  individual  parameters.  It  is  expected 
that  the  number  of  such  options  will  increase  as  the  program  is  used  for  a 
greater  number  and  variety  of  problems. 

Typiceil  options  for  non-conventional  guns  arc  those  for  gun-boosted  rockets, 
traveling-charge  guns,  and  light -gas  guns  of  the  adiabatic  compressor  type. 
Examples  of  options  for  varied  treatment  of  individual  parameters  are  those  for 
adjusted  burning  rates  (previously  mentioned),  inhibited  propellant  surfaces, 
delayed  propellant  ignition,  variable  time-intervals,  constant  resistive  pres¬ 
sure,  and  resistive  pressure  as  a  function  of  base  pressure. 
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DISCUSSION 


No  attempt  has  been  made  here  to  preeent  a  new  and  different  Interior 
ballistic  theory.  The  objective  was  to  devise  a  convenient,  flexible  scheme 
for  performing  the  tedious  numericel  calculations  required  to  obtain  detailed 
interior  ballistic  trajectories.  The  selection  of  a  program  for  high-speed 
digital  computers  has  made  it  possible  to  eliminate  most  of  the  Bintpllfications 
of  theory  required  to  facilitate  mathematical  solutions  by  other  methods. 

The  theory  presented  as  the  basis  for  the  computer  routine  is  well-known 
and  has  only  been  modified  to  account  for  coi^slte  charges.  There  are  several 
problems  present  in  all  interior  ballistic  calculations  and  these  also  prove 
troublesome  here.  For  example,  useful  propellant  burning  rates  are  not  gen¬ 
erally  available.  It  is  known  that  burning  rates  obtained  from  experimental 
firings  In  closed  chambers  are  usually  low.  The  results  obtained  from  limited 

(u) 

gun-firings  by  the  authors  '  indicate  gun  burning  rates  may  be  twice  closed 
chamber  burning  rates  under  certain  conditions.  As  previously  mentioned,  op¬ 
tional  methods  of  adjusting  closed  chamber  burning  rates  have  been  provided 
for  in  this  program.  One  such  approach  is  to  consider  the  burning  rate  to  be 
a  function  of  the  projectile  velocity  (and  possibly  a  function  of  the  projec¬ 
tile  displacement)  in  addition  to  its  known  dependence  on  pressure.  This 
method  results  in  the  use  of  closed  chamber  burning  rates  when  the  gun  chamber 
is  practically  a  closed  chamber  (v  and  x  are  effectively  zero).  When  the  pro¬ 
jectile  is  moving  at  hi^er  velocities  and  is  further  down  tube,  reasonable 
Increases  in  burning  rates  are  obtained  and  used.  Other  equally  important  dif¬ 
ficulties  are  associated  with  the  determination  of  reasonable  values  for  resis¬ 
tive  pressure  and  shot-start  pressure. 

Considerable  versatility  has  been  built  into  the  program.  Instead  of 
stopping  the  computation  at  the  end  of  Phase  III,  a  new  problem  can  be  auto¬ 
matically  read  into  the  computer  and  solved.  This  multiple-case  feature  can 
be  employed  to  advantage  for  any  number  of  additional  problems  during  a  single 
computer  run. 

Typical  Interior  ballistic  problems  were  used  to  compare  results  obtained 
from  this  computer  routine  with  results  from  other  Interior  ballistic  schemes, 
(15),  (l4),  and  (15).  The  agreement  was  generally  very  good  when  the  other 
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DISCUSSION 


No  atteiupt  has  been,  made  here  to  present  a  new  and  different  interior 
ballistic  theory.  The  objective  was  to  devise  a  convenient,  flexible  scheme 
for  performing  the  tedious  numerical  calculations  required  to  obtain  detailed 
interior  ballistic  trajectories.  The  selection  of  a  program  for  hi^-speed 
digital  computers  has  made  it  possible  to  eliminate  most  of  the  slinplificatlons 
of  theory  req.uired  to  facilitate  mathematical  solutions  by  other  methods. 

The  theory  presented  as  the  basis  for  the  computer  routine  is  well-known 
and  has  only  been  modified  to  account  for  composite  charges.  There  are  several 
problems  present  in  all  interior  ballistic  calculations  and  these  also  prove 
troublesome  here.  For  example,  useful  propellant  burning  rates  are  not  gen¬ 
erally  available.  It  is  known  that  burning  rates  obtained  from  experimental 
firings  in  closed  chambers  are  usually  low.  The  results  obtained,  from  limited 
gun-firings  by  the  authors  indicate  gun  burning  rates  may  be  twice  closed 
chamber  burning  rates  mder  certain  conditions.  As  previously  mentioned,  op¬ 
tional  methods  of  adjusting  closed  chamber  burning  rates  have  been  provided 
for  in  this  program.  One  such  approach  is  to  consider  the  burning  rate  to  be 
a  function  of  the  projectile  velocity  (and  possibly  a  function  of  the  projec¬ 
tile  displacement)  in  addition  to  its  laiown  dependence  on  pressure.  This 
method  results  in  the  use  of  closed  chamber  burning  rates  when  the  gun  chamber 
is  practically  a  closed  chamber  (v  and  x  are  effectively  zero) .  When  the  pro¬ 
jectile  is  moving  at  hi^er  velocities  and  is  further  down  tube,  reasonable 
increases  in  blaming  rates  are  obtained  and  used.  Other  equally  important  dif¬ 
ficulties  are  associated  with  the  determination  of  reasonable  values  for  resis¬ 
tive  pressure  and  shot -start  pressure. 

Considerable  versatility  has  been  built  into  the  program.  Instead  of 
stopping  the  conputation  at  the  end  of  Phase  III,  a  new  problem  can  be  auto¬ 
matically  read  into  the  computer  and  solved.  This  multiple-case  feature  can 
be  enployed  to  advantage  for  any  number  of  additional  problems  during  a  single 
conputer  run. 

Typical  interior  ballistic  problems  were  used  to  compare  results  obtained 
from  this  conputer  routine  with  results  from  other  Interior  ballistic  schemes. 
(13),  (l4),  and  (15).  The  agreement  was  generally  very  good  when  the  other 
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schemes  were  fairly  sophietleated.  In  addition,  detailed  interior  ballistic 
trajectories  are  produced  in  considerably  less  time  than  it  takes  to  calculate 
maximum  pressui-e  and  muzzle  velocity  by  other  systems.  A  typical  computer 
solution  for  a  conventional  gun  takes  only  10  seconds  if  magnetic  tape  output 
is  used  with  the  BRLESC. 

Results  from  con^uter  elmulations  have  also  been  compared  to  experimental 
data  obtained  from  well- instrumented  gun  firings.  To  demonstrate  the  adequacy 
of  the  computer  routine,  data  from  a  typical  105mm  Howitzer  firing  were  pro¬ 
cessed  by  the  method  described  in  Reference  (3.1.).  In  Appendix  D  these  experi¬ 
mental  results  are  compared  with  the  predicted  results  obtained  from  a  siiou- 
lation  of  this  firing. 
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FORM  nraCTlpN  EQUATIONS 

Oeoinetrlcal  Eiiuatloias 

1.  Initial  Volume  of  a  Propellant  Grain 


V  -  X 


(A-1) 


where;  V  =  volume  of  an  unbumed  propellant  grain,  in." 


=  outside  diameter  of  grain,  in. 

N^  =  number  of  perforations,  dimensionless 
d^  =  diameter  of  perforation,  in. 

=  length  of  grain,  in. 

2.  Volume  of  a  Partially  Burned  Propellant  Grain 

•  ==1^  “  f  [<°1  “  '^i  "  ^'i  ^^1  -  ""i^ 


2) 


where;  z.  =  mass-fraction  of  1  th  propellant  burned  at  a 
given  time,  dimensionless 

u.  =  two  times  the  distance  each  surface  has  receded 
at  a  given  time,  in. 

5.  Initial  Surface  Area  of  a  Propellant  Grain 

2  „  .  2 


r  D  "  -  N, d.^  I 

\  *  ”A)  ^  J 


(A-3) 
2 


Where;  S  =  surface  area  of  an  unbumed  propellant  grain,  in. 

D. 


4.  Surface  Area  of  a  Partially  Burned  Propellant  Grain 

Si-“  {[(”i  -  "i)  *  "A  *  “li  ][ 


(A-4) 


where:  S.  a  surface  area  of  partially  bjpned  i  th  propellant 
grain  at  a  given  time,  in. 
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Eauations  for  Nevtoa-Raphson  Method*  for  Finding  Approximate  Values  of 
the  Real  Roots  of  a  Mumerlcal  Eauatlon 

1.  Rearrange  Equation  (A-2)  to  set  f(uj^)  =  0; 

-  [  L^(N^-l)  -  S(YN^cl^)  ] 

-  [2L^(D^  +  N^d^)  + 

2.  Differentiate  Equation  (A-5)  with  respect  to  u^: 

d  [f(ui)]  f  p 

f .  (u  )  =  — i - Li  =  I  i  3(N  -l)u  2 

du^  ^  ^ 

■  2  ■  2  <»i  *  ”A)]  «j 

-  ■  "a">]} 

5.  The  value  of  the  root  of  Equation  (A-2)  Is  then: 

fCu^) 

Vi "  '^i  “  f'u^) 


(A-5) 


(A-6) 


(A-7) 


where:  =>  the  improved  value  of  the  root,  where  the 

first  estimate  of  the  root  is 


Procedure 


1.  For  each  propellant,  determine  by  integration  of  Equation  (29).  In  the 
initial  calculation  of  each  z^.  Equation  (A-3)  is  used  to  compute  each 
^i  ^®i  “  '^i  “  subsequent  calculations  of  each  z^.  Equation 

(A-4)  is  used  with  u^^  determined  as  described  below. 


*  See  Reference  (16)  for  a  discussion  of  this  method. 
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2.  The  z^'s  obtained  from  Eeiuation  (29)  are  used  to  contpute  the  u^’s  from 
£q.uatlon  (A-7)  and  then  the  new  S^'s  are  coniputed  from  Eq.uatlon  (A-4). 

In  the  initial  calculation  of  u^,  the  first  estimate  of  its  value  is  zero* 

Equation  (A-?)  is  used  to  calculate  the  improved  value,  With 

as  the  estimate,  Equation  (A-7)  is  used  again  to  calculate  a  further-improved 
value,  This  procedure  is  continued  until  the  improvement  la  lese  than 

lO'^  inch. 


APPENDIX  B 


Congputatlon  Routine 

1,  Generalized  Plow  Diagram 

2.  POBAET  Listing 


ItgTERIOR  BALLISTICS  PROGRAM  FDR  GUNS 


Interior  BallleticB  Program  for  Guns 
FORABT  LISTING 


PSOB-  1  .6^AMa — wUl.T  l^liB/llg  BUM  MAU-.-lBT-TflS - - - 

feL0G(Prtl-PB2n)  ACl-AC5)Kl-K7)YJ-Y7)Ti-Ti5)lin-lj5)Sl«S5)Rl-PS) 

hL0G(B(:i-HcSiGCl-cni»)0(:i-0f:P)F:ci-Fi;S)aci-p,r5)  1CI-1C5)  jci-jeiiiFCl-FC5i 

fiofs 

00U4 

BLnC(  G1  -(j7)Cl-C51F1-f  5 1  8A l-fi A 5 )  CO  VI -C  GY  5 )  T  0 1  -  TO 5 )  RH03  -HM0  5)  BETl -BET  5) 
BLon(nc7i-i)n7'))ANi-ANi  1.4) 

ODDS 

1  A 

B1 

CONlALPl-ALPblDl-US) tDPl'0PS)Ll-L5)NPl -NPS)XCl-XCpn)MCl-MC5) 

iHtr.Hi  A.BPAniANj  )!?)» 

U006 

7 

HEA0-FuRmAT{O1)-(WP)Ah)Vo)AP>PE  jPEl  )PPhAX)* 

PEaO-F (iRmat  (01  )  -  (Cl  )H  )GA  I  )  TO  I  )X 

8 

1008 

Bl.l 

BF.AD-PORMAT(01)-(DT)NI)KV)kX)D)EVP))(  SET(  J=0U! 

PPAO-PORmAT (01 )-( XC1.J)H«1  , J) *  COUNT (PO ) 1N( J 1  SOTO  1  HI . 1 ) * 

9 

iri 

B2 

fNTEP(  |NTEaEk)Ni)lV)KSET(  J=l))*  I N  T  (  NCP  t -PAN  )  * 

HEAn-PoRNAKOt  )-(Cl,.))Fl,.))GAl,0)COV1  .JlTOl.JlPHOI  .  J)« 

0011 

12 

>'tAU-n)Ri',(,m)i  )-(gtri,  jiALPi.  j)Ui  ,j)ripi  ,j)Li .  JiNPi.j)* 

COUNT ( .N) INI j)aOTO(H?)»SFT(J=n)* 

1012 

0  013 

B2.1 

f;NTEP(  A.PlINCH)  ANDl  )%  KNT  FP '(  A  .  PUNCH  )  ANA  V  )  1 )  * 

14 

f.NrFt.(A.PuNCH)ANP)l  )% 

15 

PIInCh-FQOhAK  0^)-<l  X  WPIXM)  VU  )  APIOFI.  1  PF  1  PPMAX  »<&>«( 

ENTEA  (  A  .PukCH)  ANBO)  1  )*  EN  1  f  R  (  A  .  Pll  NCH  1  AM  7  )  p  )  * 

16 

17 

B3 

PUNCh-fOWHAT)OJ)-<1>IC1  IFIlGAIirOlXlSNITERAA* 

Pll.'fCH-FQRM4T(04  )-<l  XCl  .  J)  FI.  J)  GA1  ,  JinnVl.  J)  T01  .  J)HH01,  J)  <A>» 

1  017 

18 

CguN  1  (  ,N)  IN(  J)l!0' 0(BA)  *  SFT(J  =  ii)*  E  NT  Fk  (  A  .  PUNCH  )  A  n89  1  1 )  * 

E  N  r  Ek ( A . PUNCH ) AN3  3 ) 1 )  * 

19 

20 

B3.1 

pi)NCk-T  OHHATlObl-O  >lB£f  1,  J)  ALPl.JlDI  .  J)DP1  .  J)L  1.  J)NP1  ,  J)<A>* 
COilNl(,NltN(J)OOTO(B3.1)»  RET(  )»[))»  FNTER(A.PUNCh)ANH9)l)* 

21' 

2? 

B3.2 

tNI£k(A,PUNCH)AN41Ji;)* 

PUNCH-FOkMAT(OH)-<]>(XCl. JlPHl, J)<A>» 

23 

24 

COUNt (POl 1M( J)fiOIO( B3.2>»  SET( J»0 )»  t  NTER ( A . PUNCH ) AN89 ) 1 ) * 

EN  rEk(  A  .PUNCH  )  ANi>7)  2)  *  8ET(SHPsBlfl.l),)Psn)STuCKs818.5)* 

25 

26 

B<t 

PUNCk-F0RHAT(n7)-<1>(DT)NnXV)KX)rVP)n)<A>* 

PC  1  bE  I  «C  I /(  SAI-D*  ACI  sbC  I /TO  I»CC  I  sF  I  ♦(;  I/TOI*  EVMsEVP4l2* 

27 

0  03  5 

84.1 

PCI  .  J  =  F1,  J»Cl,UXtftAl,  J-1 1»  AC1,.JsHC1  ,.)/T01 .  J«  CCl  ,  JbF  1 .  J4  Cl  .  J/ 
COMTTOl . J* 

0036 

0037 

nci  .o*Cl..j/RH01.J*  EC1.J  =  C1.J*C0V1 , j»  FCl . J»NP1. j-l* 
GC1.U»ll,J<NPl,J-l)-2(Dl.J+NPl,J*DPl.J)* 

0038 

0039 

HC1  ,  jsptl.l.  Jtni .  J*NP1,  J»l)Pf,  J)  +  (01 ,  J*»?-NPl,  JaDPJ  ,  J*42)* 

ICl. J=l 1. J(Dl,U*»2-NPl, J*DPl. Ja*2)« 

0040 

0041 

B5 

JC1  ..1  =  3,1  418AIC1.  J/4XC0UNT  (  .N)  INIJjGOtOIBA.DKSET  ( J«U  )* 
CTsOATPl=n* 

0042 

0043 

“IT.l 

rpi»ci .  j+GAi.j+TPi*  (;t=ci,j*ct*  comnk .n i  in(  j)goto(B5.i >* 

r,APBlP1/C'x5£T(J=n)»(;AFsQAP/(GAP-1)»FP  =  CT/wrA 

0044  ■' 
0045 

EP1 =1 +FP/0EL»TP1  =  1/ (GAP-1 )»fP 2=1/ ( (2*TP1  +3 1 /OEL  +  ( BaTP 1  +  2 ) /EP > * 
MCTD=Wk*CT/DEL*TP4=U* AGW=AP*38ft .4/wP* 

0046 

0047 

f:R2i>E)(P<0AFALoG(  i-TP2)r*  TPibExp(  i .  54i.0G(b)  )*  TPS»'EXP(S,i7i,BLaG 

0048 

CONTID))* 

004? 

B9.2 

Tp3=P  XpT;  83T?*roBTCT) )  % 

TPbbCI, J*T01. J+TP4*  COunT( .N)  IM J)G0Tn(BS.2)»  SET(JbO)* 

0050 

0051 

HCl.«(  .384l2*fPl(  XM+VO/AP)  (TPA/Ct-298))/((  1  +0,6*TP?/Tf>S) 

Tol)!" 

C0NTEVm*»2)» 

0  0  53 

BA 

nLEAR(7)N0S.AT(Kl)*  Cleah(  7)1yOS.At(y4  )*  tLPAR  (  7 1  I»jOS  ,  AT  ( 01 )  * 

0054 

KHsPKbpbhbx I «Al PB IN IHRbV*  (BUT* 

55 

B6.1 

CLEAR(5)N0S.AT(U1)*  XuiiTsn*  PRLSTsn* 

Y3,J«1.1*  COUNT(5)IN(J)UOTO(B6.1)«  SFT(Jen)* 

0056 

0057 

BS.2 
87. 1 

y3.J>0x  CnUNtl  .N)  IN(J)QC)tO(a6.2)x  SETljanix 

IF_IMT(N»l)Q()Tn(B7.!>)l|  SET  ( SHS-OR.X  .  1  )* 

o658 

0059 

IF-1nT( NB2)aOTO<B7.iH(  SeT ( BwabORS  .  1 ) »  ' 

IF-lNT(NB3)aOTO(B7.7)»  SET < SW5«DR3 . 1 ) * 

6 

OU^I 

IF-tRt( NB4)a0tO(B7.8)«  S£t(SWA>DR3,1 )X  G&tO(B3kx 

0062 

) 

I- 


45 


B7.e 

SET(£W6»814)*QOTO(Be)» 

—liHI 

BS 

TP1 an* 

0067  1 

B8.1 

TPl»rJCl,J*rPl*  COuNT(,N)IN<J>QOTO(Be.l)ll  SET(J  =  n)» 

PTsCI»f  I/O/n-TPD*  SF.T{SW1eH15)SWfl  =  R14.5)«  PMiVsPT* 

0068  1 
0069  ' 

B9 

ENTEH<R.K.G.)DT)2.N)P9)V1)K1)01)»  GnTDt.SWDX 

TF  ( vs>  =  i)RnT0<R9.i  )»set(Swii»bi  oijsn  )*r.oT0«nPi  )* 

0070  ; 
0U7i  ^ 

B9.1 

BIO 

Y3  =  li(K3s(J«S1  sO*Rl  =  ri«u3sU*SET(Swll  =  BlU  )  jeo  )!(OOTO(DRi  .  1  )K 
t  F  ( Y«>  =  i  )r.oTn(Ri  (;  .1 )  *sET(skiisbii  ).jei  )*GOTn(r)Bi  )* 

7?  1 
0  0  73  ; 

BlO.l 

Y4  =  likK4an»d2=0*P2«>0*(i4  =  l<*SErjSwll«81.1  )  J»1  l*eOTO(  ,SW3  )* 

74 

Bll 

Y5sl».Khsii«S3=0»«v'Su*l'5s0*SET(SWll!iB1  ?)  Jb2)*G0T0(  .Sw4)* 

IF(Y6>  =  l)r,()TQ(Hl2.1)i‘SET(Smi  =  8l3)Je3)»G0Tn(nR1)* 

76  . 
0077  ^ 

B12.1 

Bli 

Y6  =  l».K6a0*S4  =  r,*P4-|,‘|((;.4-u»Sf  KSWllnBl  s  )  Ja3  )  XGOTO  ( ,  Sn5  )  * 

IF  (Y7>  =  l)GflT  (1(814.  •)»SEl<SmiBRl4)J»A)*Gf)T0(nfil)» 

Y7  =  lYKV«0»S5  =  0*R5  =  ,.»O7s0*SET<SWll!:Bj4)ja4)*Gr)T0(,Sw6)* 

SET  ( jsf:  )*TP1*0» 

6o 

0001 

614.1 

TPI  *liCl  ,  J<  1  .Y3.  J)  *EC1 .  J*Y3,  J*TPl*f;OUNT(  ,N)  IM<  J)a0T0(B14 . 11* 

VC  =  V(,*Ap»XI-Tp3)IS£7  (  JaO)»TplB8CI« 

0082  ; 
0083  i 

Bl<t.2 

TPI  bHCi  4J*r3.  J  +  TPlxCOyNl  t  .NMM(J)GnT0(8l4.P)« 

sn  ( j*o  )*tp2=ac[* 

0  084 
0088 

814.3 

TPpaACl ,J*r3.J»TP2*COUNT( ,N)IN(J)a0T0(R14,3)* 

SEr(J  =  (i)*rf  HP=(TPi-ALP)/TP?*TP1  e(^(31* 

0086 

0087 

814.4 

TPl aCCI , J*y3# J»1HIkOUUNI (.NIIMIJ) 8070(814,4)* 

SFT(  jefi  )*PTaTEMP*TPl/VC*  G0T0(.SP8>* 

D08S  1 
0u69  ! 

1 

814. S 
ORl 

ENTEP( P.K.GDI* 

R1  .  JsBl  T)  .  .l♦eYP(  AI.Pl  ,  j*I.OG(PT  »  )*tJ0Bll1  ..I»HISFC1  ,  J*H2aRCl .  J« 

0090  1 
0092  ' 

i 

H3aMCl ,  J*rt4aici  ,  J»M 5s JCl .  J  ( 1- Y 3  ,  J )  *H6  at.  1  ,  J*H 7  aDl ,  J* 
hBsDPI  .JtRpBNPI  .J»hli  aDCl,  ,  J*G070(SAmP)» 

0093  : 
0094  i 

OR3 

Rl. JbRI <J*KVAK?*KX*Xl» 

K3.  JsSl»,(ARl<  J/UCl/ J* 

1094  ' 

95 

I 

0R3.1 

DR4 

GOTOJ.SwlU* 

PBBPT/FplXpRhaHB/Ep?* 

1095  : 
0096  1 

KlaAUW<P8-PHl*K2«Y1  »XIbY2« 

IF(XI<XC?n)Q0T0(UR5)»  PH*PH20*  GOTOCnPOl* 

0097  1 
0099 

i 

DRb 

DR9 

ENTEH(U,D. 1N»XI )PR)XCi)PHi)2U)3)1)3 )» 

DELX*X  I-XLSTKSa«l»RR  +  PHL3T* 

OlUO 

0098 

' 

lNTPW»( OELxiSljMl )/2*INTPH»XUST*XI*PRLSTbPR» 

ALPBl  M(;Ti)*K2**2/77p,6)*AP*INrPR+H0L*K?A4  2*  SOTO(  B1  4. 5)* 

0099 

OlUO 

1 

815 

815.1 

IF(Pl<H£)ftnT0(Bl5.1 )»  S£T(Swa3DH4)SWJ  sB16)* 

PRRbPT*  PB=Pr« 

0101 

0102 

1 

816 

lF(Yl?n)GOTO(8l7)*If ( Y3A»l)ANDe Y4>bi ) aND< Y5>sl) AND(Y6>il) 
C0NTAND( Y7>»1 )GOTO(B16.1)«GUTO(B17)» 

103 

1103 

'  j 

1 

816.1 

817 

SET(STuCKbH22)* 

XF  =  Xl/l2*V  =  Y1./l2*AFaKl/t2»SET(  jao  )»STBn* 

2103 

0104 

817.1 

B17.2 

DC71i  J  =  K3. J*C1, J»CnuNT( .N) IN( j)QOTO(Bl 7.1)*SET< J=0)8 
STsSl,J*SrxCO|lNT (.N)IN<J)a0T0(B17,2)»SET(j8n)» 

0105 

0106 

i 

B17.3 

■[  m  jn||  Jl  ■■  j  ■  IMl  M  M  U  JM  JJtf  MMj 

■tm 

B9i 

GOTO(,SWP)* 

tlY  f  £«(  A  .PUNCO  )  AND  1  )* 

loe 

109 

1 

1 

eie.2 

feNTtlR(  A.p(jNCH)AN89)l)*  TMbTbIoOO* 
PUNCM-fORRAT(Ofll»<l><TM>XI)PBR)PT)PB> V)AF  )<A>* 

110 

111 

L_ 

B18.3 

PUNCH-F0PMAT109l-<l>(TM)XI)Xf ) TE MP ) VC ) PR > ST Xfl > * 

PUNCH-FORMAT ( 01 0 ) -<1 >«  rrt)X! )Y3.J Inez J , J ) R1 . J )S 1 » j )<4> * 

112 

1112 

1 

1 

COUNT!, N)IN(J)G0T0<818.3»»  SET(J«0>* 

2112 

w  i 

I 

I 

i 

s 

i 

i 


I 


r;niiM(i8.wrPitM(.iPuliiTninio.«u  SrT(<;LiP  =  Hin  i  i.iP  =  ninnTn(  .siiini^nf 

^  o 

’ 

Bie.4 

SF:T(SWP»tH8.2)XGOTO(  .STOCK) 

4112 

B18.5 

TaT*nTX 

113 

lE(Xl>XM)G()T0tB21)»XILS1cXI»VLST  =  V«PPI  RT  =  PB* 

0114 

ENTER(R,K.(il)* 

0115 

B2] 

VmaX!:((XM-XILST)«V-VIST)/1XI-XILST))'V|ST* 

0116 

PRmAX  =  (<XM-XII  ST)(P8-PBLST)/tXt-XII  ST ) )+PRl  ST* 

0117 

! 

TPHAX  =  TPMAX*10Ii0* 

1117 

ENTEP<A.PllNCM)ANh9)l)*  K.NTE  H  (  A  .  Pi  INCH  )  A  N81  )  1  )  * 

118 

r - 

PIJNCH-E0RHAT(011)-<1>< VMAX) PMAX) XpHAX)! PMAX) PBMAX)<6>* 

119 

fiOTOlNEJRN)* 

1119 

, 

B32 

t.NTeR(A,PUNCH»ANB9)l)*£NfER(A,PUNCH)AN<57)l)*ENT£R(  A.PUNCH)AN81)1)* 

ISO 

VMAX:fl*PKMAX=n*TPMAX=TPMAX*lfin(l< 

121 

P||HCH-F0hMAT(O11)-<1>( VMAX)PMAX)XPMAX)TPMAX)PBMAX)<0>* 

122 

' 

NEWRN 

GOTOIHI )* 

123 

GAM2 

Tl=MV-U0*T2aHfl-UO 

EFl  .1 

EU=.  78  54(110A*3  +  M]-nO**2»M7»UOAH3  +  H/l)-H8* 

Bik^H 

■ 

EE  2 

FPU» .7R54(3*U0»*2*H1-2*U0*H2-H3)* 

nnrfi 

1 

EF2.1 

U01  =liU-Fii/EPl)* 

RiCfH 

EF3 

IF-ARSt  (U01-UO)<«‘.(i0001)GOTO(FF4)*|lOeliniX 

0131 

FF3.1 

OOlOiaAM?)!* 

013P 

EF4 

Sl«3 . 1416( ( m6-UO) (h7-U0  +  H9( MOtiJO) )*.5*T1aa2-.5*H9 

cnNT*T;>**?)*si .  j«si*Min/Hii»iii ,  jsijo*GOTn(nR;u* 

01 

FORM(1ii-)u)1-7) 

1 

02 

FORM!  1  2-2-9  )).l  )1?.i.1i,)1.;s)3.2)r2-l-fl)3-1  )  1  2-6-H  )  3-5  )  1  2-A-P-9  ) 

1 

03 

f  0RM( l?-?-9) 12-7-9) 3- 2) 12-1-7) 3-13) 12-4-A.25) 

134 

04 

FORM! 12-2-9) 12-7-9)3-?) 12-1 -7)3-4) 1 2-P-7) <-2) 12-4-6) 3-4)  1  ?-n-8-Sn) 

135 

05 

FORM! 12-9)3-1 ) 3  2-6) 3- 4) 12-6)3-4)17-6)3-4) J  2-1-7) 3-5)1  2-1-3-23) 

136 

06 

F0RM( 3-20)12-3-8)3-11 ) 12- 5-7-32) 

137 

07 

FORM (12-7) 3-5) 12-1-3) 3-4) 12-9)3-1) 12-9) 3-5)) 2-4-6) 3-5)12-1-7-17) 

mam 

OR 

F0RM(i2-2-H) 3-1 ) IP-3-8) 3-?) 12-6-1 n) 12-6-1 „) 12-6-1 n)l?- 5-9) 3-1) 

C0NTi2-4-1(i-9) 

ISSm 

ov 

F0HM(12-2-fi)3-l  )  12-3-8 ) 3-2) 12-2-8) 1-2)1 2-X-8) 3-2) 12-6-10) IP-4- 7)3- 3) 

Bcni 

1 

CONTlP-S-9-in) 

noa 

_ an. 

F  OHM  (  12-2-8  )  3-  1  )  Jp-3-8)  3-?)  12-1  -  7  )  :i-.1  )  1  2-4-9  )3-1  )  1  2-3-8  )  3-9)  1P-S-9-PD) 

mBm 

8  105  MM  MOwITZEH  BU  7fi5 

IPRO.l.  HT.  HAMREL  CMAMbFR  HORE  ABEA  P-K  SS  PRESS  MAX 

1  Ml  PROPELLANT 

1  CMAHEE  rORCE  fiAMMA  COvULuME  f'LAME  TEMP  DENSITY 

1  BETA  alpha  ij.U,  ShAiN  DJA.  PERe  SR.IENgTM  NO.  PERF. 

_1 _ _ _ RESISTANCE _ 

1  PHOJ.  TRAVEL  pressure 

_1 _ P-ISCELI  ANEOUS _ _ 

1  01  NO.  PROP.  nV  KX  EST.  MIJ7.  VEL.  DIAMETER 

_1 _ P  greater  THAN  DEStRED  MAX  PRESSURE _ 

1HU2ZLE  vEL,  MAX,  PHESSuRE  X  AT  PMAX  T  AT  PM4X  MljZ  PRESSyRE 

_1 _ 

I  projectile  StDPPEI; 

33. _ SIj _ iJL3_! _ 13.77 _ 46nn. _ 3.084  SOI 

.042V  1152000.  1.25  20u0, 

1-03 _ li _ g _ 0 _ 4.1.34  ISOO.  _ 

4800. 


33. 

,042V 

1-03 

.00 

_ aUL. 

.30 
,35 


APPENDIX  C 

Input  and  Output  Data 

1.  Input  Data 

2. '  Output  Data 

5«  Sample  of  Output  Format 
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1.  INPUT  DATA 


Units 

Program 

Symbol 

Gun  Constants 

Weight  of  Projectile 

lb 

WP 

Length  of  Gun  Tube 

in. 

XM 

Empty  Volume  of  Chamber 

in.5 

VO 

Cross-sectional  Area  of  Bore 

in.^ 

AP 

Shot-Start  Pressure 

pel 

PE 

Pidduck-Kent  Constant 

dimensionless 

DEL 

Resistive  Pressures 

pel 

PR1,J 

Travel  of  Projectile  Corresponding  to  each 
of  20  Resistive  Pressures 

in. 

XC1,J 

Diameter  of  Bore 

in. 

D 

Propellant  Physical  Constants 

Weights  of  Propellants 

lb 

C1,J 

Weight  of  Igniter 

lb 

Cl 

Densities  of  Propellants 

Ib/in.^ 

RHOljJ 

Outside  Diameter  of  Propellant  Grains 

in. 

D1,J 

D-i.ameter  of  Propellant  Perforations 

in. 

DP1,J 

Length  of  Propellant  Grains 

in. 

L1,J 

Number  of  Perforations  per  Grain 

dimensionless 

HP1,J 

Number  of  Propellants 

dimensionless 

N1 

Propellant  Thermodynamic  Cuuulauts’^ 

Forces  of  Propellants 

in. -Ib/lb 

F1,J 

Force  of  Igniter 

in.-lb/lb 

PI 

*  See  Reference  (1?)  for  these  data. 
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f 


Units 

Program 

Ssrmbol 

Hatloa  of  Specific  Heats  of  Propellants 

dimensionless 

GA1,J 

f 

Ratio  of  Specific  Heats  of  Igniter 

dimensionless 

OAI 

Covolume B  of  Propellants 

in.Vlt 

cbvi,j 

Adiabatic  Flame  Temperatures  of  Propellants 

°K 

1‘01,J 

Adiabatic  Flame  Ten^ierature  of  Igniter 

°K 

TOI 

Burning  Rate  Coefficients 

in.  1 

■  psi“ 

BET1,J 

Burning  Rate  Exponents  (a*s) 

dimensionless 

ALP1,J 

Burning  Rate  Velocity  Coefficient 

In. 

sec  in. /sec 

KV 

Burning  Rate  Displacement  Coefficient 

in. 

sec-ln. 

KX 

Miscellaneous  Constants 

Time  Interval 

sec 

DT 

Estimated  Muzzle  Velocity 

ft  /sec 

EVP 

Maximum  Allowable  Breech  Pressure 

psi 

PPMAX 
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2.  OUTPUT  DATA 


Identifying  Data 

The  ooniplete  llefc  of  Input  data  Is  printed  out  to  permanently  identify 
the  coniputatlon. 

Program 


Units 

Symbol 

Trajectoiy  Data 

Time 

mllllsec 

OM 

Travel  of  Projectile 

in. 

XI 

Travel  of  Projectile 

ft 

XF 

Breech  Pressure 

psi 

PKi 

Space-mean  Pressure 

psi 

PT 

Base  Pressure 

psi 

PB 

Velocity  of  Projectile 

ft/sec 

V 

Acceleration  of  Projectile 

ft/seo^ 

AP 

Temperature  of  Propellant  Gas 

0 

K 

TEMP 

Volume  behind  Projectile  available 
for  Propellant  Gas 

in.  5 

VC 

Resistive  Pressure 

psi 

PR 

Total  Surface  Area  of  Propellants 

in.^ 

ST 

Mass-fractions  of  Propellants  Burned 

dimensionless 

Y3,J 

Mass  B;imlng  Rates  of  Propellants 

Ib/sec 

DCZ1,J 

Linear  Burning  Rates  of  Propellants 

in. /sec 

R1,J 

Surface  Areas  of  Propellants 

la.® 

S1,J 

Unite 

Program 

Symbol 

Sunttnary  Data 

Muzzle  Velocity 

ft/sec 

VMAX 

Maximum  Breech  Preesure 

pel 

PMAX 

Travel  at  Maximum  Breech  Preesure 

In. 

XPMAX 

Time  at  Maximum  Breech  Preseure 

sec 

TPMAX 

Muzzle  Pressure  (Base  of  Projectile) 

psl 

PBMAX 
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3;s.Q0U00  ai.uQOOO  JS^.OOuOP  13«77000  3,n2Pon  46QO. 

________ 

CHARGE  rnRCE  GAMMA  COvO 


2,13560  3670150  ,  1.264(1  31.080  2433.  .  056700 


iI3PlHiRCIwfaHHBHilXXIiI>l 


BETA  ALPHA  b.D.  grain  UIA.  PEHF  GR. LENGTH  NO.  PERF. 

■  I)0Q5(i79  .84^(7  .  0476  ^1»4 


2 


50U00. 


!  4.0U0 

2800,  1 

'  _ _ _ _ 

26011.  1 

DT 

OOOlO 


NO.  PROP 


KV 

OOOOUU 


miscellaneous _ 

KX  EST,  MUZ.'vFL.  diameter 
.  0000000 _ isnn. _ 4.1340 


.0066  2U.907  .216  2364.96 


•  8  u  rj  0 

.000 

1531.26 

1531.26 

1531.26 

.3000 

.000 

.  0 1)  3  0 

2274,43 

.0193  22.765  .242  1656.74 


1.10(10  .000  .0000  2322.79  103.49^ 

lOQi _ .0328  31,903  _  ,340 


1.1000  .000  .0138  45.729  ,340 

I  *6ee  list  of  Output  Data  for  Rrogram  SymbolB  and  Units. 
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CD  CD  c;  a>  o  o\ch  i>j 


1.2U00 

1  ii  A  1) 

.  u  u  u 

2366.15 

1.03.3fl3  .0 

4032.79 

1.2uU0 

.  fi  n  n 

ib.SSg 

.379  1655.46 

1.2JU0 

.  ('no 

,  01.62 

51.055 

,i7i>  23/7.32 

1 . 3  U  n  0 

■  U  (1 0 

2919.97 

2919.9; 

291 9,97 

.00 

.  0001) 

1  . 3  j  u  n 

.  iiou 

.nun  i> 

^.346 . 1)9 

1  n  < .  95S 

.  n 

40.35. 09 

,  1 . 3  i)  (I  0 

.  U  0  (I 

,  il445 

39.549 

.422 

1654.39 

i.Auuj  .(ion  .i.ifta _ _  ■»??  g.'<6fi.7o 


i.^jon  .1.00  _ aUJ _ .oonn 


i.4nno 

.  noij 

.  oooo 

2355.59 

103.112 

.0 

4037.66 

1 .41111  0 

•  If  0  0 

..i5l4 

43.9)9 

.469 

1653.19 

1.4000 

.000 

.  021  7 

63.344 

.469 

2384,46 

.  1 .  S  U  U  u 

.non 

3729.29 

3728.29 

3728.28 

.00 

.  000  0 

l.SuOO 

.000 

.  noou 

2364 , DU 

102.953 

,0 

4040 , 50 

1.5000 

.000 

.0590 

49.690 

.  520 

1651.87 

1 

1 .5000 

.000 

.  11250 

70.407 

.  520 

2388.62 

1.6UU0 

.000 

4202.09 

4202.09 

420?. 09 

.00 

,  0  000 

1.6000 

.000 

.  0000 

2371.43 

102.777 

.0 

4043.64 

_ 1.6000 

.000 

.06/4 

53.998 

.576 

1650.41 

1.6000 

.000 

.  0296 

79.156 

,  576 

2393.23 

1.7000 

.000 

4  726  .'66 

4729.69 

4728.68 

.00 

.  0  00  0 

1.7000 

.000 

.OOOO 

237/. 98 

1 02. 582 

.0 

4047.11 

1.7000 

.000 

.0767 

59.574 

.637 

1648.79 

1.7000 

.000 

.  0  326 

86.655 

.537 

2398.32 

1.9000 

,000 

5384.73 

5312,49 

5169.1  D 

1.49 

8990.2031 

1.9000 

.000 

.  0000 

2363.69 

102.382 

4500.0 

4050.93 

I.BOUO 

.000 

.  0670 

65.753 

.7  04 

1647.00 

i.aouo 

.00  0 

.0371 

95.972 

.704 

2403.93 

1.9000 

.  0  04 

6040.79 

5959.75 

5798.89 

2.77 

17452.175 

1.9'JOQ 

.  1)04 

.  0003 

2398.72 

10?. 162 

4500. 0 

4055.14 

1.9000 

.  004 

.  0983 

72.457 

.777 

1645.03 

1.9000 

.004 

.0420 

106.166 

.777 

2410.11 

2.0000 

.  1'09 

6765.99 

6675. 12 

6494.96 

4.94 

26804.618 

SToooo 

1  (.09 

TocOT 

2393.04 

iOl  .^46 

4500.0 

4059, >6  “ 

2 . 0  0  0  0 

,009 

.1108 

79.730 

.956 

1642.85 

2.0000 

.  1)09 

.  0474 

117.?96 

,8  56 

2416.91 

2. 1000 

.016 

7565.33 

7463.83 

7262. 38 

8,o9 

37115.875 

2.1000 

.016 

.0018 

2396.65 

101  ,739 

4500.0 

4064.83 

2. 1000 

.  Ul6 

.  igoS 

677594 

.942 

1640*46 

2.1000 

.016 

.  0534 

129.452 

,942 

2424.37 

2.2000 

.  029 

8444.01 

8330,73 

8105.88 

12.30 

48449.250 

2.2000 

.  o26 

TooT^ 

TIW75? 

1 0 1 .856 

4500.0 

ioToTSe 

2.2000 

.029 

.1396 

96.069 

1.035 

1637,83 

2, 2  fob 

ToSl 

142.686 

1.035 

24S2.5S 

<p>)[o«  o«  Cmoi  o«  o«  m  mow  c*«i'  ck'  cmI  r\; pV  fohc  (c| 


3.8UU0  .i'Sfii?  4715,458  3,11V  0690  ,^6 


4.SUU0 

5.861 

.46/3 

559.430 

4.6000 

4 . 6 II  U  0 

6 . 558 

6. 558 

3S7#,2 .85 
.  5465 

334119.87 

2161,52 

4 . 6 11  U  0 

4 , 6  U II  0 

6.558 

6. 55B 

.997/ 

.4941 

268.442 

561.529 

4.7000 

V  .305 

32683. 8ii 

32245. 30 

4.7000 

4.7000 

7  4  A  0  S 

7  4^0^3 

,6 lift  7 

1  ,  0  0  11 U 

-tfl  3^. 

.OUC 

4 . 7ut]0 

7.305 

.  52iifc 

556.704 

4 . 8  0  U  0 

4  .  SoflO 

8  4  Id  il 
rt  4  i  U  li 

3l 54ft .35 
.6  75u 

31123.11 

2110.52 

4.8000 

4 . SOUfl 

8.1  OU 
8,100 

1  .  OOOll 
.  5469 

.>100 

544.755 

4.9000 

8.953 

3u3H2.32 

29974.7(1 

4. 9000 
4.9000 

8.94;^ 

.7452 

1  .  IIOIIU 

2085.45 

.  illi  0 

4. 9000 

8,942 

.  5724 

531.861 

8.0000 

8.0000 

9,8?9 

9, 029 

2''Pll*6.i 
.81  VI 

28819.72 

2060.93 

5.0000 

5. 0000 

0.829 

9,829 

1  .  0  0  0  (1 
.  5973 

.  OllO 
510.326 

5.  10 110 

1  II.  7  59 

2‘i:!5?.88 

276/6.52 

S.loOO 

5. louo 

10.759 

10.759 

*8966 

1  .  fl  0  M  U 

/  0  3  /  .  1  5 
•  1)  0  0 

5.1000 

10.759 

.8215 

504.441 

5.2000 

5.2000 

11  .  729 
11.729 

26920. 3o 
.9774 

26559,5  3 
2014.27 

T.  2u0  0 
5.2u0a 

11.729 

11  . 729 

1  .  iiOoO 
.6451 

.  ono 

49U . 446 

5.3000 

12.739 

25820.59 

254/4. 1/ 

5.3000 

5. 3  lion 

12. 7 
12.739 

1.  0616 

1  .  II  0 II 0 

I99<i.  j  b 
•  lUlO 

5.3000 

12.  /39 

.66/9 

4'76 .475~ 

TTTiTlTo  1  •  VTb 

S .  4  u  0  0  1  ■< .  ?  9 


a47AU.?4  2443a. 04 

t  ■  l'l8H  _ T970.H6 


3.453 


2R5«.42 


334in.fl4  505.74  023868.9(1 


166.904 

658,5 

4362.40 

3,4  36 

1480 .40 

3.436 

2882.00 

31375. on 

642,69 

414921 ,6r' 

17^.  i!7h 
.  II  f)  1.’ 

374.8 

.00 

2904.82 

3.38  0 

29(14 .82 

3  0  283.1  (. 
185.947 

682.55 
29  11.7 

402787,92 

2926.94 

,  0  (1 11 
3.28  3 

,O0 

2926.94 

1 

29165.68 

V  2  (1 . 9  5 

387359,05 

196. 4  91) 

.non 

4^6*6 

.  no 

2948,17 

3.182 

2948.1/ 

i 

28041.87 

207.686 

767.55 

883.7 

368342. 9n 

2968.51 

.fiO'i 

3 .  r.e  r. 

,00 

2968,51 

26929. 53 

792.17 

348393,83 

219.504 

.non 

1000.0 

.00 

2967.98 

2.977 

2987,98 

2684?. ?9 
231.916 

625.32 

1000.0 

333785.59  ; 

3006.61 

.  oon 

2.877 

.  OO 
3006.61 

24786.62 

8  67. 1)8 

319601.38 

24^ .  BO-* 

,  dn d 

1000. 0 

.  liO 

^024.4? 

3024.42 

23768.73 

258.417 

887.50 

1  0  0  U  .  0 

305924 . 79 

3041.46 

■>/]  l\. .6"  gatl37!)  ■  n.'< 

iiiiin.n  iOVj.ij 


5 . 8  u  11 0 

18.316 

1 .  S.'b? 

1  ^^0,^  ,  7H 

1  7  ,  b  p  0 

1  J  i,  il  .  0 

3102.53 

6 . 8 II  (1  (1 

1  P  .  .M  b 

1  .  hilfU 

.001) 

.  <M)  1* 

.  no 

ii .  8  u  11 0 

18..U  6 

.  n-yj 

4lO . 6i  / 

BDi 

31/'?  •  63 

6 .9uOO 

2-)13''.lS 

19864 . 02 

1  9.327 .89 

l!l?1  ..3  7 

>.9iinn 


1. .  ft  /<7  1  1  M76 , 4S 


1  ■!  1  1  fi  ii  n ,  0 _ 


ini^  .?i 


3ii(s.?i 


6.0IJU0  ?('.v^6  l./.'OS  3859. 8:< 

6.01100  3  66 _ 1  .  ■.i|,n. _ .  IHIQ 


>M9.774  iniHi.O  31.29.33 

I. n  I) _ _ 


6.2uQn  _  1/H9y.36  1766?,  l7l7i>.B7  1||89.34  >17341.96 

6.2000  23. .<2/  1.9439  1  826.63  38.^.6S7  IDiill.O  3164.01 

6 . 2n  0  0 _ 23.0?/ _ Lt.*JU!i9 _ _ _ .  00  0  _  .  li  p  t. _ .  0  Q _ 

6.2U00  23.32/  .rtOiie  "366.384"  2.ri43  3l.5‘i.()l 


S  .  0  46 


6 , 9  U  0  0  .4  3  .  II 4  6 


7. QUOD  34.519 
.OilOn  34.519 


7.OU00  34.519 

7.0000  34.519 


302.6u6 


322IS.79 


13467.67  13206.98  12928.37  1235.34  160271.72 

.a/66  1726.  57  533.110.3  1000.0  3235.76 


3. 0000  .000  .OOP  .00 


19500.65 


7.300 


/.3uaO  39.048 
fl 


7.3000  .39.048 


7.4000  40.59? 

7^4  00  0  __  40 *59 2 


7.4000  40.592 


l.UOOO  .000  .non  ,00 


115«1.1«  11425.80  11117.4?  1293.77  135939.51 

27  1673,65 


1 . 0  0 II 0 


7.5000  4?. 153  110.54.26  l«935.52 


7.5000  42.153  3.5127  1657.49 

7.5000  42.153  l.ojLuO  .000 


/.50UO  42,153  l.UOOO 


636,645  1000.0 

.0  00 


.000  .00 


9529.70 


.00 


.61100  43.726  10619.56  10477.09  10194.31  1319.24  123536.40 

7.6UQG  43.728  3.6440  1641.86 _ 655.238  lOOii 


7,6000  43.725  l.OQIlU  .  00  0  ,  0  00  ,  0  0 

Z._6o0,!1. 43  ._72|.  . .  1  .OfiJlQ.^ _ <0  00 _ .0011 _ .00 

7.7000  4  5.319  _ Ijj  1 84 .55  1iin47.9l  97/6.7?  1331,1  0  117925.56 

7.7000  45.319  3.7766  1.626.75  680.036  1000.0  .00 

7.7000  45,319 


7. 7000  45,319  l.UOOO  .000  .nOO 


U.gUQO 

6.5000 

6.5000 

6.5000 


5«.<t9<) 
5H.  49<1 
56.494 
56.494 


/5?e.3?  74?b.34 

4.6745  15gl.73 

1.0000  .000 

l.iiOliO  .000 


7??4.93  1409.39  83639.325 

660. Bn7  lono.o _ 

.000 

.non  .00 


6,6000  6n.]9u 


/273.04  71/5,46 


6.6000  6:), 1  90  i.onn 


8,7000  61.696  I.UOltO 

6  ■  7  u  u  u  61.696  I  •  i)  u  I'l  V 


6.6000  63.611 
8T5UOO  63.411 
6,8000  63.611 


6805.98  6/14.6/  6533.44  1432.71  74348.324 


63.611  l.UOIlO 


1  .  >100  6 


r31.067  1000  .0 

,noo  .00 


.000  ,00 


|3 


8 . 9  U  D  n 
a  .  9  :h)  0 

8 . 9  ij  n  fi 
8,9i)00 


V  .  0  0  U  0 
9.Quno 
9  .  U  (I  (i  G 
9 . 1)  U  U  H 

9 . 1  'i  0  IJ 

9.1000 

9.1000 

9.1G0U 


9 .2u0U 
9.2000 
9.2000 
9.2u0n 

9 . 3  0  U  0 
9 .  j  0  0 
9 . 3  0  0  II 
9 . 3  0  0  0 


9.4000 
9 . 4  u  0  0 
9.4000 
9.4000 


9  . »)  II I)  0 
9  .'Soon 
9  .  S  0  0 II 
9.SOOO 


9.6000 

9.6uOO 

9.6000 

9.6000 


9 . 7  u  0  0 
TmTo" 

9.70UQ 

9 . 7  u  0  0 

'5T  8  0  0  0 

9.fl000 

9.8000 

9.8000 


ZZLE.  VF 
1493.5 


6  8  .  3  4 

6890.41 

6801  .  97 

6326. Kn 

1  ,  90 

71 867.887 

68.  .*:i4 

5  .  '1  4  b 

1  477. 9il 

9  84.741 

1  0  0  0  .  0 

.  00 

A  S  t  ^  iS  H 

1  •  'i/l  •(  0 

*  II  ll  1) 

.  li  n  I) 

.  no 

68,  .<34 

1  .  ..1 0 II 0 

.  II  n  u 

,  iinii 

,  no 

A  /  «  |i6  6 

oShS.  7,3 

IJ  (1 .  y  A 

6 1  3  0 . 0  '2 

1446,83 

68927.888 

67.166 

5 . 8889 

146/. 69 

9  78. 83.3 

innn.n 

.00 

67. -66 

1  .  'U'lU 

.  imo 

,  I'Ol! 

,  no 

67 . 166 

1  .  ■.Ifliio 

.  0 II  0 

.non 

,  00 

68 . 8|i  / 

619  3 .2-/' 

8943.3  6 

1483,50 

66419. 1)46 

eehse^b 

8. 7.3  39 

148/. 78 

1 0  n? , 44  o 

1 0  0  0 . 0 

.  liU 

68.60/ 

1  .  i.  OiiO 

.  non 

.  I'fJ  ll 

.00 

6H  .  n  / 

1 ,  'innu 

.0011 

4  lip  f’ 

,00 

7  0 , 6  8  8 

6  0  ■■6.23 

8928.03 

14s>9.94 

640,32,7.34 

7n  .  ■’88 

8.fl/9(. 

144b.n7 

1020.487 

1 0  n  II ,  ll 

.  on 

7n . 888 

1  .  '1  1)  J  u 

.  DUO 

.  un  n 

.  uU 

7f.  * 

1  .  0  9  0 

.  0  II  (I 

.  f!p  n 

...zy, 1 II  _ 

8830.00 

8/81 .86 

8896.62 

1466 . 1 4 

61761 . u 28 

72.. 1111  6.0289  :438.64  lOSi.'iRI  iniiO.O  .00 

72.31n  I.JQ.iO _ .  000 _ .1.00 _ JtJJ _ 

72. . Mu  t,ii()cju  .non  .100  ,1)0 


7/1.;  8662.2/  SS«6,3l  34.'.‘'.R.<  14/2,13  S9S96.A29 

74.ii7.1  6.1/28  )429.48 _ 107  4,807  10  0  0.0 _ .00 

74 .  !■  /i  t ,  OOI'U  .  000  ,  lino  .00 

7<i .  ■  '/  3  I  ■  *  0  ;■  0  »  0  0  0 _ .  0  0 _ .  i)  0 

74. 0  4  < _ SSi.2.2/  S428. 4  4 _ 8?»-l.64  14/V.91  S7S32.803 

78.84.3  6. .3203  1420.49  1  099.1  3?  1  00  0  .  0  .  00 

7  8.043 _ LsJ/lLILO _ MIO  il  .00  I’ _ .  0  0 _ 

78.84.3  l.Jo'iO  .rmO  .lOii  .00 

^7.620  ^14 1^.5/  8277750  8138.  36  148.3. 4‘9  8SS63.3.3? 

77.620  6.4683  341  1.76 _ 1173.3  83  1  U  II  0 . 0 _ .  00 
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